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(2)  Objectives 

The  aim  of  this  proposal  is  gain  fundamental  understanding  of  the  chemical  reactivity  of 
complex  systems.  Specifically,  the  proposal  consists  of  two  experimental  projects:  (i)  the  state- 
correlation  of  product  pairs  in  chemical  reactions  of  polyatomic  molecules,  and  (ii)  the  solvation 
dynamics  in  water  clusters  and  in  aqueous  solution,  in  particular  the  ion-solvent  interactions  that 
are  fundamental  to  many  important  phenomena  in  chemistry  and  biology. 

(3)  Status  of  the  efforts 

The  crossed  molecular  beam  project  is  an  established  one  in  this  laboratory.  In  2008  we 
took  the  full  advantage  of  our  unique  capability  of  measuring  the  product  pair  correlation  and 
continued  making  significant  contributions  to  the  field  of  chemical  reaction  dynamics,  as 
evidenced  from  the  invitations  to  write  an  article  to  Proc.  Nat’l  Acad.  Sci.  USA  and  to  several 
international  conferences  as  a  plenary  speaker. 

The  solvation  dynamics  is  a  new  project  funded  mainly  by  Academia  Sinica  from  2008. 
Several  initial  tasks  were  carried  out,  including  modifications  of  an  existing  apparatus  for  water 
cluster  generation,  the  purchase  of  a  state-of-the-art  femtosecond  laser  system,  and  the 
installation  of  a  clean-room  facility  to  accommodate  the  new  laser  system  etc.  Those  time- 
consuming  preparation-phase  works  are  mostly  completed  by  the  end  of  2008. 

(4)  Abstract 

We  continued  the  fundamental  study  of  the  reactivity  of  Cl-atom  towards  methane.  Being 
a  competing  reaction  to  Cl  +  O3  and  one  of  the  major  sinks  for  CH4  (a  greenhouse  gas),  this 
reaction  plays  a  crucial  role  in  atmospheric  chemistry  and  is  highly  relevant  to  the  ozone 
production/depletion  problems.  Our  aim  is  to  understand  how  different  forms  of  reagent  energy 
(translation  and  vibration)  affect  the  reaction  rate  and  detail  dynamics.  Following  our  earlier 
studies  on  Cl  +  CHD3(vi=1)  and  Cl  +  CH4(v3=l),  we  now  extend  to  Cl  +  CH2D2(vi=1  and  V6=l). 
Here,  vi=l  and  V6=l  indicate  one-quantum  excitation  of  the  CH2-symmetric  stretching  and  CH2- 
antisymmetric  stretching  mode  of  the  CH2D2  reagent,  respectively.  These  two  modes  are  nearly 
degenerate,  but  with  very  different  vibrational  motions.  Because  the  oscillation  strengths  of  both 
IR-transitions  are  relatively  weak,  the  experimental  challenge  is  to  pump  up  as  much  CH2D2 
molecules  as  possible  into  the  desired  vibration  mode.  To  this  end,  we  have  designed  and 
fabricated  a  multipass  ring  reflector  (Rev.  Sci.  Instrum.  79,  033105  (2008))  that  effectively 
increases  the  infrared  laser  power  by  10  folds,  thus  making  the  experiments  feasible.  Preliminary 


analysis  of  the  data  reveals  some  very  exciting  results,  which  could  drastically  shape  our 
thinking  about  mode-  and  bond- selective  chemistry  in  the  future. 

In  the  other  project  aimed  to  understand  the  solvation  dynamics  of  aqueous  solution,  some 
preparatory  progresses  have  also  been  made.  This  is  a  new  project  funded  by  a  very  competitive 
5-year  Academia  Sinica  Investigator  Award,  starting  from  2008.  The  major  parts  of  a  state-of- 
the-art  femtosecond  laser  system  have  been  purchased  and  just  arrived.  In  order  to  accommodate 
the  new  laser  system,  the  whole  machine  was  shut  down  and  a  clear-room  unit  installed.  We  are 
now  in  the  process  of  re-installing  the  machine  and  integrating  it  with  the  new  laser  system. 
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We  report  a  comprehensive  study  of  the  quantum-state  correlation 
property  of  product  pairs  from  reactions  of  chlorine  atoms  with 
both  the  ground-state  and  the  CH  stretch-excited  CHD3.  In  light  of 
available  ab  initio  theoretical  results,  this  set  of  experimental  data 
provides  a  conceptual  framework  to  visualize  the  energy-flow 
pattern  along  the  reaction  path,  to  classify  the  activity  of  different 
vibrational  modes  in  a  reactive  encounter,  to  gain  deeper  insight 
into  the  concept  of  vibrational  adiabaticity,  and  to  elucidate  the 
intermode  coupling  in  the  transition-state  region.  This  exploratory 
approach  not  only  opens  up  an  avenue  to  understand  polyatomic 
reaction  dynamics,  even  for  motions  at  the  molecular  level  in  the 
fleeting  transition-state  region,  but  it  also  leads  to  a  generalization 
of  Polanyi's  rules  to  reactions  involving  a  polyatomic  molecule. 

mode-specific  chemistry  |  Polanyi's  rules  |  product  pair  correlation 

One  of  the  central  concepts  in  understanding  chemical 
reactivity  is  the  idea  of  the  transition  state  (1).  When  two 
reactants  collide  to  form  reaction  products,  old  bonds  are  broken 
and  new  bonds  formed.  This  bond  breaking  and  forming  process 
must  then  occur  over  a  molecular  configuration  intermediate 
between  reactants  and  products,  and  this  intermediate  region  of 
potential  energy  landscape  is  what  chemists  refer  to  as  the 
transition  state.  The  transition  state  is  often  located  near  the  top 
of  a  potential  barrier  and  acts  as  a  bottleneck  in  a  chemical 
reaction.  The  properties  of  the  transition  state  therefore  deter¬ 
mine  the  reaction  rate,  as  well  as  many  of  the  more  detailed 
observables  from  scattering  experiments. 

Over  the  past  decades,  there  has  been  tremendous  progress 
in  experimental  characterization  of  the  structure  of  the  transi¬ 
tion  state,  notably  by  using  the  spectroscopic  probes  (2-4). 
Transition-state  spectroscopy  experiments  performed  to  date 
are  essentially  the  half-collision  type  in  which  the  transition  state 
is  directly  accessed  either  through  photodetachment  of  negative 
ion  precursor  in  a  frequency-resolved  experiment  (3)  or  by  the 
femtosecond  pump-probe,  time-resolved  approach  (4).  As  ele¬ 
gant  and  informative  as  those  experiments  are,  half-collision 
results,  in  general,  do  not  depict  a  full  picture  of  how  the 
reactants  transform  into  the  products.  One  way  to  think  of  this 
is  as  follows.  The  basic  idea  of  a  typical  half-collision  experiment 
is  to  initiate  the  reaction  at  transition  state  by  a  photoexcitation 
process.  By  virtue  of  photoabsorption,  the  total  angular  momen¬ 
tum,  that  is,  the  partial  wave  or  the  impact  parameter,  of 
the  reactive  system  is  then  well  specified  and  often  limited  to  the 
lowest  few  quantum  numbers  in  a  restricted  geometry  of  the 
Franck-Condon  region.  Consequently,  the  half-collision  results 
are  greatly  simplified  and  more  amenable  to  theoretical  tests.  In 
contrast,  a  chemical  reaction  inevitably  constitutes  the  contri¬ 
bution  from  collisions  with  a  full  range  of  impact  parameters  and 
orientations.  The  resultant  wave-interference  patterns,  arising 
from  the  coherent  sum  of  scattering  amplitudes  of  many  partial 
waves,  are  manifested  in  the  full-collision  attribute  such  as 
product  angular  distribution  (5,  6),  which  cannot  be  readily 
accounted  for  by  the  few-partial-wave,  half-collision  approach. 
On  the  horns  of  a  dilemma,  a  full-collision  experiment  usually 
deals  with  asymptotic  properties  of  the  reaction,  thereby  ren¬ 
dering  direct  probes  of  the  fleeting  transition  state  difficult. 

Here,  we  propose  an  approach  to  delineate  the  dynamical 
aspects  of  the  transition  state  in  a  full-collision  experiment  by 


tracking  the  energy  flow  along  the  reaction  path.  We  previously 
introduced  an  experimental  method  to  unfold  the  state-specific 
correlation  of  coincident  product  pairs  in  polyatomic  reactions 
(7-9).  More  recently,  we  exploited  the  product  pair-correlation 
approach  to  elucidate  mode-selective  chemistry  of  the  Cl  + 
CHD3(vi  =  1  or  v3  =  1)  ->  HC1  +  CD3(v  =  0)  reaction  (10).  In 
the  latter  study,  we  found  that,  contrary  to  the  current  percep¬ 
tion,  C-H  stretch  (vi)  excitation  is  no  more  efficient  than  an 
equivalent  amount  of  translational  energy  in  enhancing  the 
reaction  rate;  CD3  bend  (v3)  excitation  is  only  moderately  more 
efficient.  These  unexpected  results  then  raised  an  important 
question:  How  do  we  understand  the  observed  differential 
reactivity  between  polyatomic  reactant  vibration  and  translation 
from  the  perspective  of  Polanyi’s  rules  (11,  12)?  The  work 
reported  here  presents  an  all-important  complement  to  resolve 
the  apparent  paradox  by  mapping  out  the  complete  energy-flow 
pattern  through  correlating  as  many  coincidently  formed  prod¬ 
uct  pairs  as  possible  to  an  initially  prepared  reactant  state.  What 
emerged  is  a  conceptually  appealing  picture  in  which  the  coop¬ 
erative  motion  of  atoms  in  passing  through  the  transition  state 
can  be  visualized.  In  addition,  this  conceptual  framework  leads 
naturally  to  a  generalization  of  Polanyi’s  rules  to  a  reaction 
involving  polyatomic  molecules. 

What  Are  Polanyi's  Rules? 

Simply  stated,  Polanyi’s  rules  concern  how  the  barrier  location 
influences  the  energy  requirement  and  the  energy  disposal  in  a 
direct  atom  +  diatom  chemical  reaction  (11,  12).  For  an  exo¬ 
thermic  A  +  BC  reaction,  the  reaction  barrier  is  usually  located 
in  the  entrance  valley  of  the  reaction,  that  is,  an  early  barrier. 
According  to  Polanyi’s  rules,  reactant  translational  energy  is 
then  more  effective  than  vibration  to  surmount  the  barrier  to 
reaction,  thus,  accelerating  the  reaction  rate.  The  converse  will 
be  true  for  an  endothermic,  late-barrier  reaction.  By  the  prin¬ 
ciple  of  microscopic  reversibility  (1),  the  total  available  energy 
will  then  be  deposited  mostly  into  product  vibration  for  an 
early-barrier  reaction,  whereas  a  translationally  hot  product  will 
be  yielded  from  a  late-barrier  reaction.  Hence,  the  rule  eluci¬ 
dates  the  role  of  different  forms  of  energy  (vibration  versus 
translation)  in  an  elementary  chemical  reaction,  and  links  its 
intimate  relationship  to  the  underlying  feature  (the  barrier 
location)  of  the  three-atom  interaction  potential. 

Experiment  on  Product  Pair-Correlated  Images 

We  performed  the  experiment  under  single-collision  conditions 
by  using  a  crossed  molecular  beam  apparatus  (7,  10,  13,  14).  A 
discharge-generated,  pulsed  Cl  beam  (5%  Cl2  seeded  in  He  at  6 
atm  total  pressure)  was  double-skimmed  and  directed  to  cross 
with  a  pulsed  CHD3  molecular  beam  in  a  high-vacuum  chamber. 
A  tunable  infrared  (IR)  optical  parametric  oscillator/amplifier 
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Fig.  1.  Raw  images,  with  (Left)  and  without  ( Center )  IR-pumping,  of  probed 
product  states  from  the  Cl  +  CHD3  reaction  at  Ec  =  8.1  kcal/mol.  REMPI  bands 
used  to  probe  the  methyl  products  are  indicated  in  parentheses.  For  CHD2(1!) 
no  detectable  signals  can  be  observed  from  the  IR-off  image.  Also  exemplified 
for  the  CD3  states  are  the  three  IR-off  images  acquired  at  Ec  ^16  kcal/mol 
(Right),  that  is,  with  approximately  the  same  total  energy  as  the  C-H  stretch- 
excited  reactant  at  Ec  =  8.1  kcal/mol.  Because  of  the  weak  signals,  some 
backgrounds  (appearing  as  blurred  spots)  were  observed  but  discarded  in  data 
analysis.  The  ringlike  feature  can  be  ascribed,  on  energetic  grounds,  to  the 
product-state  pair  as  labeled  (see  text  for  notations). 


prepared  the  CHD3  reactant,  before  the  collision  center,  with 
one-quantum  excitation  along  the  C-H  stretching  bond  via  the 
vi  =  0  — »  1,  R(l)  transition  at  3,005.57  cm-1  (15).  After  the 
collision,  the  reaction  products,  either  CD3  or  CHD2  radicals, 
were  probed  by  (2  +  1)  resonance-enhanced  multiphoton  ion¬ 
ization  (REMPI)  spectroscopy  ^331-339  nm  depending  on  the 
REMPI  bands  (16-18),  and  a  time-sliced  ion  velocity  imaging 
technique  mapped  the  state  correlation  of  coincidently  formed 
coproducts  HC1  or  DC1  (7-9).  (Under  the  experimental  condi¬ 
tions  of  this  study,  the  estimated  scaling  (up)  factors  for  probing 
the  2i,  22,  and  23  states  of  CD3  radical  are  9.0  ±  0.5,  3.5  ±  0.5, 
and  16.4  ±  2.0,  respectively.  The  notation  of  2*  refers  to  the 
vibrational  mode  2  (the  umbrella-bend)  with  /-quantum  excita¬ 
tion.)  Pair-correlated  state  and  angular  distributions  were  then 
exploited,  in  light  of  available  ab  initio  theory  (19-22),  to  unveil 
the  microscopic  reaction  pathways.  We  further  sharpened  the 
comparison  with  the  result  of  a  ground-state  reaction  at  either 
the  same  initial  translation  energy  ( Ec )  or  the  higher  Ec  with  an 
equivalent  amount  of  total  energy  (vibration  +  translation). 

Fig.  1  presents  several  raw  images  with  the  probe  laser 
frequencies  fixed  at  the  peak  of  the  Q  branch  of  the  respective 
REMPI  bands.  The  vibronic  band  notation  in  the  figure  is  such 
that  2 1  designates  the  spectroscopic  transition  involving  the  v2 
(umbrella-bend)  mode  with  one-quantum  excitation  each  in 
both  the  electronically  ground  state  (the  subscript)  and  the 
electronically  excited  state  (the  superscript).  Methyl  product 
state  was  probed  with  the  IR-excitation  on  and  off  in  an 


alternative  manner.  Very  small  signals  from  the  C-H  stretch- 
excited  reaction  were  also  detected  for  one-quantum  excitation 
in  the  v4  (in-plane  bend)  and  the  v3  (antisymmetric  stretch) 
modes  of  the  CD3  product,  but  they  were  too  weak  to  be 
quantified.  No  other  CD3  or  CHD2  modes  showed  detectable 
signature  from  the  C-H  stretch-excited  reaction,  although  some 
of  them  indicated  very  weak  signals  from  the  ground-state 
reaction.  It  is  quite  remarkable  that,  despite  numerous  product 
states  being  energetically  accessible,  so  few  product  vibrational 
modes  are  active  in  this  reaction.  Typical  signals  for  vibrationally 
excited  methyl  products  shown  in  Fig.  1  were  only  a  few  percent 
of  that  for  the  ground-state  product.  In  addition  to  the  less 
favorable  REMPI-detection  sensitivity  when  probing  the  vibra¬ 
tionally  excited  products  than  the  0q  origin  band  (16),  the  loss  in 
signals  from  the  C-H  stretch-excited  reaction  also  arises  from 
the  fact  that  only  ^20%  of  reactants  were  excited  by  the  IR  laser 
(10,  13,  14). 

The  energetics  of  the  reaction  are  well  defined:  the  reaction 
endothermicities  are  1.73  and  1.94  kcal/mol  for  the  H  and  D 
atom  abstraction  channels,  respectively.  Ec  was  kept  ***8.1  kcal / 
mol  and  the  initial  ro-vibration  excitation  of  CHD3(vi  =  1  ,j  = 
2)  adds  another  8.63  kcal/mol  to  the  total  energy.  (The  K 
quantum  number  of  CHD3  was  not  resolved  in  this  study.)  By 
conservation  of  energy  and  momentum,  the  maximum  velocities 
of  the  coproduct  HC1  (or  DC1),  recoiling  from  the  state-selected 
CD3  (or  CHD2),  in  different  vibration  states  were  calculated  and 
identified  as  the  ringlike  features  on  images  in  Fig.  1.  The  clear 
separation  between  the  rings  indicates  unequivocally  the  low 
rotational  excitation  of  the  HC1  (or  DC1)  coproduct.  The  relative 
intensity  of  the  ring  on  a  given  image  reflects  the  probability  for 
the  coincident  formation  of  the  corresponding  HC1  or  DC1  state. 
The  intensity  around  each  ring  then  gives  an  immediate  impres¬ 
sion  about  the  preferred  scattering  direction  of  the  product  state. 
[Superimposed  on  each  image  in  Fig.  1  is  a  red  arrow,  pointing 
to  the  0°-angle  that  is  defined  as  the  direction  (in  the  center- 
of-mass  frame)  of  the  initial  CHD3  beam.]  Inspection  of  the 
image  reveals  rich  variations  not  only  among  different  methyl 
product  states,  but  also  for  a  given  state  under  the  three  different 
experimental  conditions. 

Such  variations  are  better  appreciated,  after  data  analysis,  in 
terms  of  the  pair-correlated  vibration  branching  and  angular 
distribution  (7,  8).  As  indicated  in  Fig.  1,  the  inner  ring  on  each 
of  the  three  CD3  IR-on  images  (Left,  top  three)  constitutes  two 
nearly  degenerate  components.  (Recall  that  the  vibrational 
energy  of  stretch-excited  reactants  is  8.63  kcal/mol  and  the 
formation  of  HCl(v'  =  1)  requires  at  least  8.24  kcal/mol.)  A 
forward  ringlike  feature  that  is  absent  on  the  IR-off  image  (Fig. 
1  Center)  corresponds  to  the  concomitantly  formed  HCl(v'  =  1) 
from  the  stretch-excited  reaction.  In  the  side-  and  back- 
scattering  directions,  however,  the  signal  of  this  product  pair 
overlaps  with  the  contribution  of  the  HCl(v'  =  0)  pair  from  the 
ground-state  reactants  that  are  unpumped  by  the  IR  laser.  To 
disentangle  the  pair-correlated  angular  distribution  of  the 
stretch-excited  reaction  from  the  IR-on  image,  the  fraction  of 
C-H  stretch-excited  reactants  (typically  ^20%)  in  the  CHD3 
beam  was  first  determined  by  the  threshold  method  (13).  By 
scaling  down  the  IR-off  angular  distribution  by  0.20  to  account 
for  the  unpumped  ground  state  CHD3  and  subtracting  it  from 
the  IR-on  data,  the  genuine  distribution  from  the  stretch-excited 
reaction  was  then  uncovered  from  the  overlapped  ringlike 
feature.  The  results  of  such  analysis,  along  with  those  for  the 
outer  rings,  are  summarized  in  Fig.  2. 

For  the  Cl  +  CHD3(vi  =  1)  -^HCl(v')  +  CD3(vO  reaction,  all 
three  product  pairs  associated  with  HCl(v'  =  1)  display  a  similar 
angular  pattern  (Fig.  2^4):  a  sharp  forward  peak  superimposed 
on  a  nearly  isotropic  component.  The  sharpness  of  the  forward 
peak,  however,  descends  in  the  order  of  (1,  00)s  >  (1,  2|)S  >  (1, 
22)s.  Here,  the  product-state  pair  is  labeled  as  follows:  the 


12668  |  www.pnas.org/cgi/doi/10.1073/pnas.0800220105 


Yan  etal. 


0  45  90  135  180  45  90  135  180 


6cm  (de9ree) 

Fig.  2.  Summary  of  the  state-correlated  angular  distributions  at  Ec  =  8.1 
kcal/mol.  The  distributions  for  the  ground-state  reaction  at  Ec  =  16  kcal/mol 
shift  more  forward  as  anticipated  (1 0, 24, 25),  thus  not  shown  here.  To  contrast 
the  variation  in  shapes,  the  angular  distributions  in  each  panel  are  plotted 
with  approximately  the  same  peak  heights. 

numbers  in  the  parentheses  denote  (from  left  to  right)  the 
quanta  of  vibrational  excitation  in  HC1  and  the  mode  (vi)  in  CD3 
products,  respectively;  the  inner  subscript  specifies  the  quantum 
of  CD3  mode  and  the  outer  subscript  indicates  the  ground  (g)  or 
the  stretch-excited  (s)  reactant  states.  Angular  distributions  for 
the  product  pairs  associated  with  HCl(v'  =  0)  from  the  stretch- 
excited  reaction,  Fig.  2B ,  display  rather  different  patterns:  Both 
(0,  00)s  and  (0,  2i)s  pairs  show  predominantly  backward-sideways 
distributions,  whereas  the  (0,  22)s  pair  indicates  a  significant 
forward  preference. 

For  the  ground-state  reaction  (Fig.  2 C),  the  angular  distribu¬ 
tion  of  the  (0,  Oo)g  pair  displays  a  characteristic  sideways  peak 
accompanied  by  a  sharp  cutoff  against  forward  scattering, 
indicative  of  a  direct  reaction  mechanism  governed  by  large 
impact-parameter  collisions,  that  is,  peripheral  dynamics  (23- 
25).  The  distributions  for  (0,  2i)g  and  (0,  22)g  also  exhibit  the 
sharp  forward  cutoff,  albeit  more  backscattered,  suggesting  a 
direct  rebound  mechanism  with  significantly  more  contributions 
from  smaller  impact-parameter  collisions  than  the  (0,  Oo)g  pair. 
This  trend  corroborates  well  with  the  chemical  intuition  that  for 
a  collinear  Cl-H-C  transition  state,  the  smaller  impact- 
parameter  collisions  will  preferentially  lead  to  the  umbrella- 
excited  CD3  products.  Comparing  Fig.  2  B  and  C,  the  formation 
of  some  scattered  products  near  the  0°-angle  for  the  (0,  Oo)s,  (0, 

21) s,  and  (0,  22)s  pairs  (Fig.  2 B)  is  particularly  noteworthy  (see 
below).  Aside  from  these  forward-scattered  features,  it  is  in¬ 
triguing  to  note  that  the  global  shapes  of  the  angular  distribu¬ 
tions  for  (0,  Oo)s  and  (0,  2i)s  resemble  those  for  (0,  2i)g  and  (0, 

22) g,  whereas  the  distributions  for  (0,  22)s  and  (0,  Oo)g  seem  alike 
in  appearance. 

As  to  the  DC1  +  CHD2  isotope  channel  (Fig.  2D),  the 
observed  angular  distributions  for  all  product  pairs  from  either 
the  ground-state  or  the  CH  stretch-excited  reactant  are  virtually 
identical.  The  dominance  of  side-  and  back-scattered  products  is 
reminiscent  of  typical  direct  abstraction  reaction  governed  by 
rebound  mechanism  (1). 

Visualizing  the  Cooperative  Atomic  Motions  While  a  Chemical 
Reaction  Is  Taking  Place 

The  Conceptual  Framework.  To  shed  more  light  on  the  dynamics 
underlying  the  intricate  angular  pattern,  we  examined  the  rela¬ 
tive  branching  ratio  of  the  product-state  pair.  In  deriving  the 


Cl  +  CHD3 


S2  (amuv2bohr)  Si  ( amu1l2bohr ) 


Fig.  3.  Schematic  representation  of  vibrationally  adiabatic  potential  energy 
curves  along  the  reaction  coordinate  5.  (Right,  left-hand  side)  The  HCI  +  CD3 
(DCI  +  CHD2)  isotope  channel.  For  clarity,  only  those  states  relevant  to  this 
study  are  shown.  The  shaded  areas  (near  S  =  0  and  -0.5  amu1/2bohr)  denote 
the  regions  of  strong  curvature  and  Coriolis  couplings,  where  vibrationally 
nonadiabatictransitions  occur.  Forthe  HCI  +  CD3  isotope  channel,  the  relative 
cross-sections  of  different  product-state  pairs  under  three  different  experi¬ 
mental  conditions  are  normalized  and  represented  by  the  colored  bars;  those 
forthe  DCI  +  CHD2  product  channel  are  normalized  independently.  The  color 
codes  are:  green  and  blue,  ground-state  reactions  at  Ec  =  8.1  and  16  kcal/mol, 
respectively;  red,  the  stretch-excited  reaction  at  8.1  kcal/mol.  The  estimated 
uncertainties  associated  with  each  number  are  ±10%,  ±15%,  and  ±15%  for 
the  21f  22,  and  23  pairs,  respectively. 

branching  ratio  of  each  pair,  we  normalized  its  flux  to  the  (0,  Oo)g  % 
pair  from  the  ground-state  reaction  at  Ec  =  8.1  kcal/mol.  Both  i 
the  fraction  of  stretch-excited  reactants  (^20%  of  total)  and  the 
different  detection  sensitivity  when  probing  the  excited  CD3 
products  (16)  were  taken  into  account  in  the  HCI  +  CD3 
channel.  Relative  sensitivities  of  detecting  CHD2(v  =  0  and  vi  = 

1)  are  yet  to  be  calibrated;  thus,  their  normalizations  are  just 
based  on  signal  strengths  and  are  independent  of  the  HCI  +  CD3 
channel.  The  final  results  are  summarized  in  Fig.  3,  along  with 
the  adiabatic  correlation  of  the  relevant  vibrational  energy  levels 
leading  to  both  isotopic  product  channels. 

By  using  the  reaction  path  Hamiltonian  approach  (26),  pre¬ 
vious  ab  initio  calculations  of  isotopically  analogous  reactions 
mapped  out  the  minimum  energy  path  and  the  evolution  of  the 
vibrational  frequencies  of  various  modes  along  the  reaction  path 
(19-22).  By  adding  the  theoretically  calculated  vibration  fre¬ 
quencies  (with  isotope  corrections)  to  the  minimum  energy  path, 
we  connected  the  vibrational  energy  levels  according  to  their 
symmetries,  employing  the  rule  that  energy  levels  for  vibration 
of  the  same  symmetry  do  not  cross  (26). 

Theory  predicted  that  as  the  Cl  atom  approaches  the  H  end 
of  CHD3,  the  chemical  interaction  induces  a  static  curvature 
coupling  (i.e.,  coupling  of  a  vibrational  mode  to  the  reaction 
coordinate  induced  by  the  curvature  of  the  reaction  path) 
between  the  C-H  stretching  (vi)  motion  and  the  reaction 
coordinate,  resulting  in  a  significant  decrease  of  its  frequency  in 
the  transition-state  region  (19-22)  and  allowing  energy  flow 
between  this  mode  and  the  reaction  coordinate.  Similar  behavior 
was  found  for  the  CD3  umbrella  mode  (v3),  yet  other  modes  show 
little  variation  in  frequencies.  Theoretical  calculations  further 
predicted  that  these  two  active  vibrations  (vi  and  v3)  significantly 
couple  to  each  other  via  Coriolis  interactions  (21),  that  is,  the 
intermode  mixings  induced  by  the  twisting  of  the  two  transverse 
vibrations  about  the  curved  reaction  path  as  the  reaction  pro¬ 
ceeds.  Both  curvature  and  Coriolis  couplings  are  particularly 
strong  near  the  shaded  regions  in  Fig.  3.  As  the  Cl  atom 
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Fig.  4.  Classification  of  the  activity  of  reactant  vibrational  modes  in  a 
chemical  reaction.  By  using  the  present  reaction  for  illustration,  an  example  of 
the  spectator  mode  is  the  C-H  stretching  excitation  in  the  Cl  +  CHD3(\/i  =  1)— » 
DCI(v'  =  0)  +  CHD2(vi  =  1)  channel,  and  of  the  adiabatic  mode  in  the  Cl  + 
CHDbCv^  =  1)  -»  HCI(i/  =  1)  +  CD3(v  =  0)  channel.  The  bend-excited  CHD3 
behaves  as  a  transitional  mode  in  the  Cl  +  CHD3(v3  and/or  ve  =  1)  -»  HCI(v'  = 
0)  +  CD3(v  =  0)  reaction. 


approaches  the  D  atom,  however,  the  CD3  symmetric  stretch  (v2) 
is  now  brought  to  couple  to  the  reaction  coordinate  (19).  Near 
the  transition  state,  it  mixes  extensively  with  the  umbrella  mode 
(v3),  the  CD3  rock  (v6),  and  possibly  the  CD3  deformation  (vs) 
mode.  In  contrast,  the  C-H  stretching  frequency  now  becomes 
invariant,  in  accord  with  the  chemical  intuition  that  the  nonre¬ 
acting  C-H  bond  behaves  as  a  spectator  when  the  Cl  atom  attacks 
the  D  end  of  CHD3. 

Fig.  3  is,  of  course,  merely  the  vibrational  correlation  diagram. 
Vibrational  motion  rarely  behaves  entirely  adiabatically  during  the 
course  of  a  chemical  reaction.  Moreover,  in  analogy  to  the  elec¬ 
tronic  Bohn-Oppenheimer  approximation  (1),  the  concept  of 
vibrational  adiabaticity  is  rooted  on  the  relative  time  scales  of  the 
vibrational  period  and  the  interaction  time.  Although  the  latter  is 
governed  by  the  “slow”  motion  of  the  two  heavy  reactants  in  the 
present  reaction,  the  adiabatic  concept  could  become  blurred 
for  the  low-frequency  modes.  Keeping  this  in  mind  and  allowing  for 
vibrational  nonadiabaticity,  Fig.  3  serves  as  our  starting  point  for 
visualizing,  at  least  in  a  qualitative  sense,  the  energy  flow  while  bond 
breaking  and  bond  formation  are  taking  place. 

To  make  the  concept  more  concrete  and  to  set  the  stage  for 
further  discussion,  Fig.  4  classifies  the  limiting  behaviors  of 
different  vibrational  modes  in  a  chemical  reaction.  Reactant 
vibration  is  called  conserved  if  it  is  retained  as  one  of  the 
vibrational  motions  of  the  product.  During  the  reaction,  it  can 
either  remain  as  a  spectator  (i.e.,  preserving  its  mode  character 
with  the  vibrational  frequency  nearly  unchanged  throughout  the 
reaction  path)  or  behave  adiabatically  (i.e.,  preserving  the  vi¬ 
brational  quantum  number  but  with  varying  frequencies  due  to 
the  static  curvature  coupling  to  the  reaction  coordinate).  In  the 
former  case,  vibrational  excitation  in  the  nonreactive  bond 
(spectator)  does  not  actively  participate  in  the  reaction;  thus,  the 
initial  excitation  in  that  bond  is  likely  to  be  retained  in  the  final 
product  vibrational  distribution.  In  the  adiabatic  case,  energy 
exchanges  between  the  vibration  and  the  motion  along  the 
reaction  coordinate  can  occur;  thus,  an  adiabatic  vibration  is  an 
active  participant  in  the  reaction.  Another  type  of  active  mode  is 
the  transitional  vibration  (27),  for  which  the  vibrational  motion 
of  the  reactant,  usually  of  a  low-frequency  mode,  does  not 
correlate  to  any  product  vibration,  rather  it  evolves  into  the 
rotation  and  translation  of  the  departing  products. 

In  general,  the  activity  of  a  reactant  vibrational  mode  will  fall 
into  at  least  one  of  the  above  types:  spectator,  adiabatic,  and 
transitional.  Although  a  spectator  (transitional)  mode  is  always 
conserved  (active),  the  adiabatic  mode  can  partake  in  both 
behaviors  in  a  reaction.  It  should  also  be  pointed  out  that  the 
classification  here  refers  to  the  limiting  cases  and  is  not  always 
unambiguous.  A  given  mode  may  change  its  activity  from  one 
type  to  the  other  along  the  reaction  path  due  to  the  curvature  and 


intermode  Coriolis  couplings.  Consequently,  the  occurrence  of 
bifurcating  reactive  fluxes  into  multiple  microscopic  reaction 
paths  may  well  be  a  norm  rather  than  an  exception. 

As  the  Cl  Atom  Abstracts  the  Stretch-Excited  H  Atom.  A  recent  study 
on  the  ground-state  Cl  +  CHD3  reaction  demonstrated  that 
nearly  98%  of  products  were  the  ground-state  pair  (0,  00)g  when 
CD3(v  =  0)  was  probed  (10);  a  similar  result  was  obtained  for 
Cl  +  CH4  (25).  Detailed  analysis  of  the  IR-off  images  (Fig.  1) 
revealed  that  the  relative  branching  ratios  of  the  (0,  2i)g  and  (0, 
22)g  pairs  to  (0,  Oo)g  at  Ec  =  8.1  kcal/mol  (Fig.  3,  green  bars)  are 
merely  0.05  and  0.01,  respectively,  which  increase  slightly  to  0.07 
(or  0.19/2.7)  and  0.05  (or  0.14/2.7)  even  at  Ec  =  16  kcal/mol  (Fig. 
3,  blue  bars).  The  ground-state  reaction  is  therefore,  by  and 
large,  vibrationally  adiabatic.  Previous  experiments  on  reactions 
with  umbrella-excited  CHD3  (10)  or  CH4  (28)  demonstrated  the 
dominance  of  the  (0,  Oo)b  product  pair;  its  angular  distribution  is 
almost  identical  to  that  from  the  ground-state  reaction  at  the 
same  Ec ,  suggesting  instead  a  predominantly  nonadiabatic  path¬ 
way.  In  other  words,  because  of  the  strong  curvature  coupling  in 
the  entrance  valley,  the  umbrella-bending  vibration  of  methane 
actually  behaves  as  a  transitional  mode  rather  than  as  an 
adiabatic  mode  that  would  have  yielded  the  (0,  2i)b  product  pair 
by  vibrational  correlation. 

For  the  C-H  stretch-excited  reactant,  the  reaction  proceeds 
initially  over  the  v\  =  1  potential  energy  surface  at  long  range. 
As  the  intermolecular  distance  decreases,  a  (avoided)  crossing  of 
the  vi  =  1  and  v3  =  3  (labeled  as  3v3  in  Fig.  3)  adiabatic  curves 
in  the  entrance  valley  could  mix  in  some  CD3  umbrella-bending 
character  into  the  C-H  stretching  motion.  Around  Si  ~  -0.5 
amu1/2bohr,  strong  curvature  couplings  take  place,  resulting  in  a 
bifurcation  of  reactive  trajectories  into  a  vibrationally  nonadia¬ 
batic  as  well  as  the  adiabatic  pathways.  As  can  be  envisioned 
from  Fig.  3,  energy  flowing  out  of  the  initially  deposited  C-H 
stretch  will  be  greatly  facilitated,  via  Coriolis  couplings,  by  the 
proximate  umbrella-excited  energy  surfaces  near  the  shaded 
regions.  Compared  with  the  branching  ratios  for  (0,  2i)g  and  (0, 
22)g  from  the  ground-state  reaction  (Fig.  3,  green  or  blue  bars), 
the  analogous  product  pairs  from  the  C-H  stretch-excited  reac¬ 
tion  (Fig.  3,  red  bars)  show  significantly  larger  ratios,  in  support 
of  this  interpretation.  In  conjunction  with  the  strong  curvature 
couplings  of  both  (stretching  and  umbrella-excited)  active  modes 
to  the  reaction  coordinate  (19-22),  a  cascading  energy  flow 
could  lead  to  a  higher  population  in  the  ground-state  pair  (0,  Oo)s 
than  the  (0,  2i)s  pairs,  which  is  exactly  shown  in  Fig.  3. 

The  above  nonadiabatic  reaction  path  is  a  direct  mechanism, 
yet  it  is  mediated  by  the  umbrella  motion  and  invokes  a  different 
reaction  profile  from  the  ground-state  reaction.  Observation  of 
different  angular  distributions  for  the  (0,  00)s  and  (0,  00)g  pairs 
(Fig.  2)  is  therefore  not  too  surprising.  As  to  the  v2-excited 
product  pairs  (0,  2i)s  and  (0,  22)s  from  the  stretch-excited 
reaction,  additional  pathways  might  also  come  into  play  (see 
below).  Nonetheless,  the  striking  contrast  in  the  general  appear¬ 
ance,  as  alluded  to  earlier,  between  the  two  sets  of  analogous 
pair-correlated  angular  distributions  from  the  ground  and 
stretch-excited  reactant  states  (Fig.  2  B  vs.  C)  is  intriguing  and 
calls  for  further  theoretical  work  for  deeper  understanding. 

Not  all  reactive  trajectories  undergo  nonadiabatic  transitions; 
those  that  remain  vibrationally  adiabatic  will  retain  their  original 
character  in  the  sense  that  the  one  quantum  of  stretching 
excitation  is  preserved  within  the  Cl-H-C  moiety  of  the  colliding 
pair.  Those  trajectories  could  be  temporarily  trapped  by  the 
dynamic  well  associated  with  the  stretch-excited  adiabatic  curve 
in  the  transition  state  region  (Fig.  3),  allowing  more  time  for 
energy  redistribution  (10,  25).  Angular  distributions  for  the  (1, 
00)s,  (1,  2i)s,  and  (1,  22)s  product  pairs  (Fig.  2 A)  are  distinct  from 
the  corresponding  pairs  with  HCl(v'  =  0)  products  (Fig.  2 B), 
showing  a  forward  peak  on  top  of  an  isotopic  component — a 
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characteristic  distribution  for  reaction  involving  a  short-lived 
complex  (6,  29-33).  We  assert  that  those  products  paired  with 
HCl(v'  =  1)  are  produced  predominantly  from  a  complex¬ 
forming  pathway,  and  conceivably  are  mediated  through  a 
resonance  state  trapped  by  the  dynamic  well  (10,  25). 

Further  intracomplex  energy  redistribution  might  occur  before 
the  resonant  complex  decays  adiabatically  to  the  HCl(v'  =  1) 
product  pairs.  In  particular,  the  proximity  of  the  (0,  22)  curve  near 
Si  ^0  amu  1/2bohr  suggests  that  it  is  a  competing  nonadiabatic  path. 
The  inverted  branching  ratios  between  (0,  22)s  and  (0,  2i)s,  0.79  vs. 
0.50  (Fig.  3),  support  this  view.  The  observation  that  the  relative 
reactive  fluxes  in  the  forward  direction  (Fig.  2 B),  which  could  be 
regarded  as  the  imprint  of  the  vibrationally  nonadiabatic  decay  of 
resonant  complexes,  decrease  progressively  in  the  order  of  (0,  22)s 
>  (0, 2i)s  >  (0,  Oo)s  also  closely  corroborates  this  scenario.  It  is  worth 
noting  that  from  the  branching  ratios  listed  in  Fig.  3,  the  excited 
HCl(v'  =  1)  pairs  from  the  C-H  stretch-excited  reaction  (Fig.  3,  red 
bars)  collectively  account  for  ^30%  of  total  reactivity  of  the  HC1 
+  CD3  isotope  channel.  Treating  it  as  a  rough  estimate  of  vibra¬ 
tional  adiabaticity,  the  overall  nonadiabatic  pathways  appear  quite 
facile. 

As  the  Cl  Atom  Attacks  the  Unexcited  D  Atoms.  For  the  other  isotope 
channel  DC1  +  CHD2,  only  three  product  pairs,  (0,  li)s,  (1,  00)s, 
and  (0,  00)s,  are  significantly  populated  from  the  C-H  stretch- 
excited  reactant.  Theory  predicts  that  the  C-H  stretching  fre¬ 
quency  hardly  changes  as  the  Cl  atom  attacks  the  D  atoms  (19), 
implying  that  the  initial  C-H  excitation  remains  localized  as  a 
spectator  in  forming  the  adiabatically  correlated  product  pair  (0, 
li)s.  And  the  shape  of  the  C-H  stretch-excited  reaction  path  for 
forming  the  (0,  li)s  pair  should  be  nearly  identical  to  that  for 
producing  the  (0,  Oo)g  pair  from  the  ground-state  reactant  (Fig. 
3).  If  the  stretch-excited  reaction  indeed  proceeds  adiabatically, 
then  the  angular  distribution  of  the  (0,  li)s  pair  should  be  similar 
to  the  (0,  Oo)g  distribution.  That  is  exactly  what  we  observed  (Fig. 
2D).  Yet,  the  measured  branching  ratios  (the  red  bars)  showed 
significant  variance  with  that  from  the  ground-state  reaction 
(either  the  green  or  blue  bars).  In  particular,  the  adiabatically 
correlated  (0,  li)s  pair  accounts  for  only  30%  of  total  reactivity 
of  this  isotope  channel  (a  value  identical  to  the  other  isotope 
channel  may  be  fortuitous),  in  significant  deviation  from  the 
adiabatic  expectation  or  the  spectator  paradigm  (34, 35)  that  the 
initial  excitation  of  the  unreactive  C-H  bond  survives  as  the  (0, 
li)s  product  pair.  Therefore,  the  initial  C-H  excitation  is  coun¬ 
terintuitively  not  a  mere  spectator  when  a  D  atom  is  abstracted. 

The  angular  distribution  for  the  dominant  (0,  00)s  pair  in  the 
D  atom  abstraction  channel  is  virtually  identical  to  the  other  two 
pairs  (Fig.  2D).  The  formation  of  the  (0,  00)s  pair,  however,  must 
involve  a  facile  nonadiabatic  pathway  to  funnel  the  energy 
initially  deposited  in  the  C-H  bond  into  the  rotational  and 
translational  motions  of  the  departing  products.  What  kind  of 
cooperative  nuclear  motions  might  partake  in  redistributing  the 
initially  localized  C-H  stretching  energy  as  the  Cl  atom  is 
abstracting  a  D  atom?  Theoretical  calculations  suggested  that 
the  CD3  symmetric  stretch  (v2)  mode  of  CHD3  is  an  active  mode 
in  this  isotope  channel  (19).  As  depicted  in  Fig.  3,  we  conjectured 
that  several  proximate  combination  modes  involving  the  v2  mode 
of  CHD3,  for  example,  v2  +  v6,  are  the  plausible  candidates  for 
nonadiabatic  transitions,  as  illustrated  by  the  dashed  circle,  in  the 
entrance  valley.  [The  rocking  vibration  v6  is  a  transitional  mode, 
which  preferentially  leads  to  the  rotational  and  translational 
motions  of  reaction  products  (19).]  The  shape  of  the  reaction 
path  for  v2  +  v6  is  uncertain  because  of  the  strong  mixings 
between  the  v6  and  v3/vs  modes  in  the  transition-state  region, 
thus,  plotted  as  a  dashed  line  in  Fig.  3.  Tentatively,  the  reaction 
starts  with  a  C-H  stretch-excited  reactant  in  the  entrance  valley. 
Near  the  circled  region,  ^30%  of  reactive  fluxes  stay  adiabatic 
and  eventually  yield  the  (0,  li)s  product  pair;  the  other  70%  of 


reactive  fluxes  make  a  nonadiabatic  transition  to  the  combina¬ 
tion-mode-excited  path.  As  the  reaction  proceeds  further,  ex¬ 
tensive  couplings  occur  around  S2  ~  -0.5  amu1/2bohr  and  energy 
can  flow  into  the  other  degrees  of  freedom.  Approximately 
two-thirds,  or  0.45/(0.45  +  0.22),  end  with  the  ground-state 
product  pair  through  nonadiabatic  processes,  and  one-third,  or 
0.22/(0.45  +  0.22),  forms  as  the  (1,  00)s  pair. 

Mode-Specific  and  Bond-Selective  Reactivity 

These  issues  are  at  the  heart  of  polyatomic  reactivity  and  have 
been  actively  pursued  both  experimentally  (10,  28,  34-37)  and 
theoretically  (19-22)  in  recent  years.  We  have  reported  a  strong 
mode  specificity  in  terms  of  pair-correlated  distributions  for  two 
different  modes  of  excitation:  the  C-H  stretch  and  umbrella 
bend  of  the  Cl  +  CHD3  HC1  +  CD3  reaction  (10).  As 
presented  above  (Fig.  3,  the  colored  bars),  sharp  contrasts  on  the 
pair-correlated  branching  ratios  in  both  isotope  channels  are  also 
noted  when  compared  with  the  C-H  stretch-excited  and  ground- 
state  reactions.  It  is  insightful  here  to  have  a  global  view  of  the 
total  reactivity,  that  is,  the  sum  of  the  pair-correlated  branching 
ratios  in  each  isotopic  channel,  under  three  different  experi¬ 
mental  conditions. 

For  the  DC1  +  CHD2  channel,  the  total  reactivity  of  the 
ground-state  reaction  at  Ec  =  8.1  kcal/mol  (Fig.  3,  green  bar),  the 
C-H  stretch-excited  reaction  at  the  same  Ec  (the  sum  of  the  red 
bars),  and  the  ground-state  reaction  at  16  kcal/mol  (the  blue 
bars)  are  1.0,  0.94,  and  2.47,  respectively.  A  total  reactivity  of 
0.94  for  the  C-H  stretch-excited  reactant  is  not  much  different 
from  the  ground-state  reactivity  of  1.0  at  the  same  Ec.  In  other 
words,  although  the  initial  one-quantum  excitation  of  the  C-H 
stretch  (a  spectator  bond  here)  exerts  enormous  effects  on 
product  state  distributions  that  deviate  from  the  spectator 
picture  because  of  the  breakdown  of  vibrational  adiabacity,  it 
does  not  alter  much  the  total  reactivity  of  the  D  atom  transfer  | 

channel.  The  latter  conclusion  seemly  conforms  to  the  spectator-  5 

bond  paradigm  that  the  vibrational  energy  in  the  nonreacting 
bond  yields  little  effect  on  the  reaction  rate  (34,  35).  Hence, 
spectator  or  not  depends  on  the  measured  quantity. 

As  to  the  HC1  +  CD3  channel,  the  total  reactivity  of  the 
ground-state  reaction  is  1.06  (the  sum  of  the  green  bars)  and  3.1 
(the  blue  bars)  at  Ec  =  8.1  and  16  kcal/mol,  respectively;  and  the 
relative  reactivity  for  the  stretch-excited  reactant  at  8.1  kcal/mol 
becomes  4.44  (the  red  bars).  Hence,  with  the  equivalent  amount 
of  total  energy,  the  stretching  vibration  is  more  effective  than 
pure  translation  energy  by  a  factor  of  ^1.4  (or  4.44/3.1)  in  this 
isotope  channel.  This  finding  seems  in  accord  with  the  expec¬ 
tation  of  Polanyi’s  rules,  as  well  as  with  the  chemical  intuition 
that  the  vibrational  energy  is  directly  deposited  into  the  bond 
(C-H)  to  be  broken.  However,  the  enhancement  factor  of  1.4  for 
a  stretch-excited  reactant  is  virtually  the  same  as  that  obtained 
for  a  bend-excited  reactant  (10).  Thus,  the  preferential  promo¬ 
tion  of  total  reactivity  by  vibration  in  the  HC1  +  CD3  channel 
does  not  appear  to  be  mode-specific. 

Moreover,  when  both  isotope  channels  are  considered,  the 
combined  enhancement  factor  by  stretch-excitation,  (4.44  +  0.94)/ 

(1.06  +  1.0)  =  2.61,  turns  out  to  be  nearly  the  same  as  the 
translational  enhancement  of  2.7,  that  is,  (3.1  +  2.47)/(1.06  +  1.0). 

In  other  words,  vibration  is  no  more  efficient  than  translation  in 
promoting  the  overall  reactivity.  These  seemingly  conflicting  views 
about  the  reactant  vibrational  effects — more  effective  in  one  iso¬ 
tope  channel  but  not  in  overall  reactivity — stem  from  the  different 
activities  of  the  C-H  stretch-excitation  partaking  in  the  two  isotopic 
product  channels:  behaving  as  an  adiabatic/transitional  mode  in  the 
H  atom  abstraction  channel  and  as  a  spectator/transitional  mode 
when  the  D  atom  is  transferred.  Consequently,  although  the  C-H 
stretching  vibration  is  more  effective  than  translation  in  promoting 
the  formation  of  HC1  +  CD3,  the  converse  is  true  for  the  DC1  + 

CHD2  channel. 
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In  terms  of  the  overall  isotopic  product  branching  ratio 
[HC1+CD3]/[DC1+CHD2],  the  stretch-excited  reactant  increases 
the  ratio  from  1.06  for  ground-state  reactant  at  Ec  =  8.1  kcal/mol 
to  4.72  (or  4.44/0.94),  which  is  to  be  compared  with  a  value  of  1.26 
(or  3.1/2.47)  for  translationally  hot  ground-state  reaction.  This 
differential  isotope  effect  between  the  stretch-excitation  and  trans¬ 
lational  energy,  or  the  preferential  cleavage  of  the  excited  C-H 
bond,  is  a  manifestation  of  bond-selective  chemistry  (36,  37). 

A  Renaissance  of  Polanyi's  Rules  for  Polyatomic  Behavior 

How  do  we  reconcile  the  present  viewpoint  with  the  venerable 
Polanyi’s  rules  (11, 12)  that  are  taught  in  textbooks  to  understand 
the  energy  disposal  and  requirement  in  a  direct  atom  +  diatom 
reaction  (1,  38)?  Transition-state  structure  in  an  exothermic 
three-atom  reaction  is,  in  general,  reactant-like  according  to 
Hammond’s  postulate  (1);  thus,  an  attractive  surface  or  an  early 
barrier  is  predicted  by  Polanyi’s  rules.  Because  the  barrier  is 
located  in  the  entrance  valley  where  the  potential  shape  trans¬ 
verse  to  the  reaction  coordinate  is  not  yet  strongly  perturbed  by 
the  intermolecular  interactions,  the  (diatomic)  reactant  vibra¬ 
tional  frequency  should  not  alter  significantly.  Hence,  from  Fig. 
4,  the  reactant  vibration  behaves  as  a  spectator  during  motions  up 
to  the  transition-state  region,  and  little  vibrational  enhancement 
in  reactivity  is  expected.  However,  for  an  endothermic  A  +  BC 
reaction,  the  structure  of  the  transition  state  will  be  product-like 
and  the  barrier  lies  late  in  the  exit  valley.  As  the  reaction 
proceeds  and  before  the  transition  state,  the  colliding  pair  must 
pass  through  the  corner-cutting  region  where  the  vibrational 
frequency  orthogonal  to  the  minimum  energy  path  changes 
because  of  the  curvature  coupling  to  the  reaction  coordinate, 
facilitating  the  energy  exchange  between  the  initial  vibration 
excitation  and  the  motion  along  the  reaction  coordinate.  There¬ 
fore,  the  vibrational  energy  of  reactant  becomes  effective  in 
promoting  reaction  rate  and  yields  translationally  hot  products. 
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Reactant  vibration  is  now  active  and  often  behaves  as  a  transi¬ 
tional  mode  in  reaction. 

As  one  might  surmise  from  the  above,  the  premise  behind  our 
approach  to  extending  Polanyi’s  rules  to  a  polyatomic  reaction  is  the 
recognition  of  different  activities  or  responses  of  multiple  vibra¬ 
tional  modes  toward  the  interacting  atom  (Fig.  4).  Although  the 
numbers  of  the  vibrational  modes  with  motions  orthogonal  to 
the  reaction  coordinate  scale  as  3N-7,  which  can  be  quite  large  as 
the  size  of  molecule  increases,  possible  types  of  their  activities  in  a 
reaction  may  remain  quite  limited.  (Taking  the  present  6-atom 
reaction  as  an  example,  the  number  of  internal  (or  vibrational) 
degrees  of  freedom  of  the  system  evolves  from  9  on  the  reactant  side 
to  12  at  the  transition  state,  and  to  7  for  the  two  molecular  products. 
During  the  course  of  chemical  reaction,  one  of  the  12  degrees  of 
freedom  near  the  transition  state  corresponds  to  the  reaction 
coordinate  and  the  other  11  modes  are  the  vibrational  motions 
orthogonal  to  the  reaction  path.)  Illustrating  in  this  report  a 
prototypical  Cl  +  CHD3  reaction  outlines  our  attempt  to  qualita¬ 
tively  comprehend  the  dominant  factors  that  govern  the  polyatomic 
reactivity.  Theoretically,  we  adopted  the  reaction  path  Hamiltonian 
approach,  along  with  previous  high-level  ab  initio  calculations  for  an 
approximate  yet  illuminating  elucidation  of  the  mode  activities  and 
the  intermode  couplings  to  account  for  vibrational  nonadiabaticity. 
Experimentally,  product  pair  correlation  was  exploited  to  disen¬ 
tangle  the  intricate  pathways  at  the  vibrationally  correlated  level 
(Fig.  3).  Putting  them  together,  the  cooperative  nuclear  motions 
during  the  course  of  a  chemical  reaction  were  qualitatively  decoded 
and  then  unveiled.  Depending  on  the  motional  response,  different 
vibration  modes  will  exhibit  different  mode-  or  bond-selective 
behaviors.  Hence,  within  the  conceptual  framework  proposed  here, 
the  “generalized”  Polanyi’s  rules  can  be  regarded  as  a  ramification 
or  the  other  side  of  the  coin  of  mode-specific  reactivity. 
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Abstract 

Using  a  photoelectron  imaging  technique,  remarkable  variations  of  photoelectron  angular 
distributions  are  demonstrated  in  the  region  of  autoionization  resonance  states  of  a  two-photon 
excited  I-atom.  It  is  shown  that  the  observed  energy  evolution  of  angular  variation  exhibits 
Fano-like  behaviour  and  can  be  parametrized,  for  an  isolated  autoionization  state,  using  the 
theoretically  proposed  scaling  relationship  for  the  width  and  the  energy  shift  with  respect  to 
the  absorption  spectrum  of  the  autoionizing  state.  In  the  case  of  overlapping  resonances, 
however,  it  is  found  that  the  scaling  relationship  fails  to  account  for  the  energy  dependences  of 
the  ^4  parameters,  though  it  remains  valid  for  the  spectra.  Possible  reasons  are  discussed. 


1.  Introduction 

Asymmetric  lineshape  is  ubiquitous  in  many  atomic  and 
molecular  spectra  [1 , 2] ,  as  well  as  in  condensed  matter  physics 
[3-7].  The  universal  character  of  this  profile  [8]  is  recognized 
as  the  manifestation  of  a  quantum  interference  phenomenon, 
as  first  pointed  out  by  U  Fano  in  interpreting  the  autoionization 
spectra  of  rare  gases  [9].  Here,  the  process  arises  from 
the  excitation  to  a  discrete  resonance  state  embedded  in  the 
ionization  continuum.  Thus,  two  competing  optical  pathways 
are  possible:  one  path  links  the  ground  state  to  an  excited 
discrete  state,  and  the  other  connects  the  ground  state  with  a 
continuum  of  energy  states.  The  interference  of  two  alternative 
pathways  gives  rise  to  typically  an  asymmetric  profile  of 
absorption  cross  sections  as  a  function  of  energies,  which  can 
be  expressed  as  the  Fano  formula  [10,  11]: 

a  =  ab  +  oa(q  +  e)2/(l  +e2).  (1) 

Here  ah  is  the  excitation  cross  section  for  transitions  into 
continuum  that  do  not  interact  with  the  superexcited  Rydberg 
state  (i.e.  the  background  absorption  far  away  from  resonance), 
oa  is  that  for  transitions  into  the  continuum  coupled  to  the 
Rydberg  state  and  s  =  (E  —  Er)/(T / 2)  is  the  reduced  energy 
with  Er  and  T  being  the  peak  position  and  the  width  of 
the  resonance,  respectively.  The  dimensionless  profile  index 
or  Fano  parameter  q  measures  quantitatively  the  degree  of 


asymmetry  of  the  resonance  profile  from  the  interference 
of  the  two  alternative  excitation  pathways  between  resonant 
excitation  and  direct  photoionization. 

To  gain  a  clear  physical  insight  into  the  Fano  formula  (1), 
figure  1  summarizes  the  basic  idea  and  the  effects  of  the  q 
value  on  the  spectral  profile.  Depicted  in  figure  1(a)  are  the 
excitation  mechanisms  and  the  interference  pathways,  where 
Ci  and  C2  are  the  two  alternative  excitation  amplitudes  and 
Ci 2  denotes  the  coupling  amplitude  between  the  discrete  and 
continuum.  In  this  scheme,  the  Fano  parameter  q  corresponds 
toCjyC2Ci2,  thus,  measures  the  strength  of  the  coupling  to  the 
continuum.  The  decisive  role  of  the  q  parameter  in  dictating 
the  lineshape  is  illustrated  for  q  <  0  in  figure  1(b),  where  the 
relationship  between  the  q  value  and  the  characteristic  features 
(the  extremes)  of  the  spectral  profile  is  explicitly  shown.  Also 
depicted  in  figure  1(b)  are  the  cases  for  a  positive  q  and  the 
two  limiting  values,  0  and  00.  As  can  be  seen,  the  symmetric 
Lorentzian  peak  is  recovered  when  q  — >  00,  i.e.  when  the 
coupling  between  the  discrete  state  and  the  continuum  (C12) 
is  vanishingly  weak.  On  the  other  hand,  when  this  coupling 
is  strong  and  the  transition  probability  to  the  discrete  state  is 
weak,  then  q  0  and  a  symmetric  ‘window  resonance’  is 
observed. 

While  the  Fano  formula  (1)  is  well  known  and  has  been 
widely  used  to  fit  and  interpret  the  asymmetric  spectral  peak, 
it  is  less  documented  for  the  other  experimental  observables. 
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Figure  1.  Schematic  representation  of  the  basic  mechanism  underlying  the  Fano  lineshape  (a),  in  which  two  indistinguishable  pathways — a 
direct  excitation  C2  and  an  indirect  one  through  Q  and  the  discrete-continuum  coupling  Ci2 — are  possible  to  reach  the  same  final 
continuum  state.  Four  different  profiles  are  shown  in  (b)  to  illustrate  the  effect  of  the  profile  index  q  on  the  spectral  lineshape.  Also 
exemplified  in  the  case  of  q  <  0  are  the  relationships  between  the  characteristic  features  and  the  q  value. 


Since  the  absorption  cross  section  (or  the  spectral  profile) 
and  the  vector-correlation  parameters  are  merely  two  sides 
of  the  same  coin — both  manifest  the  behaviour  of  ionization 
amplitudes — one  may  ask:  When  an  autoionization  peak 
shows  pronounced  asymmetric  absorption  profile  and  how 
does  the  other  concomitant  observable  such  as  photoelectron 
angular  distribution  or  the  product  ion  polarization  behaves 
across  the  autoionizing  resonance  peak?  A  recent  work  by 
Grum-Grzhimailo  et  al  [12]  addressed  exactly  such  a  question. 
They  showed  theoretically  that  in  the  region  of  an  isolated 
autoionizing  resonance  state,  the  vector  correlation  parameters 
also  exhibit  the  Fano-like  behaviour,  and  the  width  and  the 
peak  position  of  the  vector  correlation  observables  possess 
universal  scaling  relationships  to  those  encountered  in  the 
usual  photoabsorption  cross  section.  These  predictions  were 
validated  by  the  concurrent  experimental  measurements  of  the 
alignment  and  orientation  of  the  final  photoion  in  the  region 
of  the  Xe  4d<2/26p(/  =  1)  autoionizing  resonance  [12]. 

Recently,  this  laboratory  reported  a  photoionization 
dynamics  study  of  iodine  atom  using  the  photoelectron 
imaging  technique  [13].  Numerous  autoionizing  resonances 
were  identified  and  a  variety  of  photoelectron  angular 
distributions  were  observed.  During  the  course  of  that 
investigation,  we  noticed  that  the  observed  angular  distribution 
appeared  sensitively  dependent  on  the  exact  photon  energy, 

1. e.  within  1  cm-1  of  dye  laser  output.  An  in-depth 
investigation  was  then  undertaken,  which  revealed  peculiar 
yet  systematic  trends  of  photoelectron  angular  distributions 
across  autoionizing  peaks.  In  an  effort  of  trying  to  understand 
those  intriguing  results,  we  came  across  the  above  paper  by 
Grum-Grzhimailo  et  al.  As  reported  here,  their  theoretical 
framework  provides  an  insightful  means  to  elucidate  the 
underlying  physics  of  our  observations. 

2.  Experiment 

As  the  experimental  details  have  been  presented  previously 
[13,  14],  only  a  brief  description  will  be  given  here.  A  cold 
molecular  beam  was  generated  by  expanding  a  gas  mixture 
of  ~9.5%  CH3I  in  He  at  a  total  pressure  of  30  psi  via  a 
fast  opening  Even-Lavie  pulsed  valve.  The  molecular  beam, 
which  was  skimmed  and  collimated  to  ~1  mm  diameter 
along  the  ion  time-of-flight  axis,  was  intersected  at  right 


Figure  2.  The  photoelectron  spectrum  (a)  and  the  photoionization 
spectrum  of  the  fragment  I-atom  (b)  following  the  photodissociation 
of  CH3I.  The  numbers  label  the  spectral  peaks  of  the  iodine  atom 
(see  [13]). 

angle  by  a  frequency-doubled  laser  beam.  The  laser  was 
vertically  polarized  (i.e.  parallel  to  the  detector  face)  and 
focused  by  an/ =  50  cm  spherical  lens  to  the  centre  between 
the  repeller  and  extractor  electrodes.  In  this  one-colour 
experiment,  the  laser  photolyzed  CH3I  and  ionized  either 
the  parent  molecule  [14]  or  the  neutral  photofragments — 
mainly  the  I-atom  [13].  Electrons  produced  in  the  interaction 
region  were  probed  by  the  velocity-map  imaging  technique 
[15,  16].  The  non-uniformity  of  the  detection  sensitivity  of 
the  imaging  system  was  corrected,  following  the  procedure 
detailed  previously  [13].  The  anisotropy  coefficients,  and 
@4,  of  the  photoelectron  angular  distribution  were  derived  from 
fitting  the  Abel-inverted  images  [17]. 

3.  Results  and  analysis 

3.1.  Autoionization  spectral  profile 

To  set  the  stage,  figure  2  displays  two  photoionization  spectra 
by  (a)  monitoring  the  photoelectron  and  (b)  detecting  the 
iodine  cation  as  in  previous  reports  [13,14].  The  peak  numbers 
label  (some  of)  the  atomic  iodine  transitions  over  this  energy 
range.  As  discussed  previously  [13],  those  peaks  arise  from 
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Figure  3.  Four  representative  photoelectron  images  at  two-photon 
energies  near  peak  6,  with  (a)-(d)  marked  in  figure  4.  In  each  panel, 
the  left  half  displays  the  raw  image  and  the  right  half  shows  the 
Abel-inverted  one.  The  laser  polarization  axis  is  vertical  in  the 
plane.  The  arrows  point  to  the  ring  feature  arising  from  the 
autoionization  process  of  the  I-atom.  The  other  ring  structures  can 
be  assigned  to  the  vibronic  autoionization  of  the  CH3I  molecule  [14]. 

two  types  of  photoionization  processes:  the  (2+1)  resonance- 
enhanced  multiphoton  ionization  (REMPI)  process  of  the 
ground  state  I(2P3/2)  atom  and  the  two-photon  transitions 
originating  from  the  spin-orbit  excited  state  I(2Pi/2)  and 
terminating  in  superexcited  Rydberg  (or  autoionization)  states 
I**  that  converge  to  the  I+(3Pi)  ionization  threshold.  The  latter 
peaks  are  of  concern  in  this  study,  namely,  the  autoionization 
process  of 

I(2Pi/2)  +  2 hv  [I**]  -*  I+(3P2)  +  e“.  (2) 

Figure  3  shows  four  photoelectron  images  acquired  at  different 
laser  wavelengths  in  the  vicinity  of  peak  6.  In  each  panel,  the 
left  half  displays  the  raw  image  and  the  right  half  shows  the 
Abel-inverted  one  [13,  15].  On  energetic  grounds,  the  multi¬ 
ring  features  can  readily  be  assigned  to  the  photoionization 
of  either  the  parent  molecule  CH3I  [14]  or  the  I-fragment 
(indicated  by  the  arrow)  [13],  from  which  the  fraction  of  the 
latter  intensity  to  the  total  intensity  can  then  be  determined. 
To  recover  the  photoionization  spectrum  of  our  interest, 
process  (2),  from  the  total  photoelectron  spectrum  shown  in 
figure  2(a),  numerous  photoelectron  images  across  a  given 
spectral  peak  were  taken  and  the  dependence  of  the  intensity 
fraction  on  photon  energy  was  then  determined.  Figure  4 
shows  such  results  for  peak  6.  Displayed  in  the  upper  panel 
are  the  energy  dependences  of  the  total  electron  intensity 
and  the  fraction  coming  from  process  (2);  the  resultant 
autoionization  spectral  profile  is  presented  in  the  lower  panel. 
Also  shown  in  the  lower  panel  is  the  photoionization  spectrum 
of  the  same  peak  by  monitoring  the  cation  I+  (figure  2(b)). 
Clearly,  the  I+-spectrum  is  considerably  narrower,  suggesting 


that  other  processes,  such  as  multiphoton  ionization  of  CH3I 
followed  by  dissociation,  are  taking  place  [13]  and  contribute 
to  the  observed  I+  signals,  thus  distorting  the  spectral  profile. 
Note  by  passing  that  figure  3  gives  a  clear  demonstration 
of  the  power  of  the  ion  (or  electron)  velocity  map  imaging 
technique  in  unravelling  the  desired  ionization  process  from 
other  concomitantly  occurring  processes. 

3.2.  Isolated  resonance 

The  angular  distributions  from  the  autoionization  process 
(2),  indicated  by  arrows  in  figure  3,  exhibit  a  dramatic 
dependence  on  photon  energies.  For  example,  it  changes 
from  predominantly  positive  /3  values  in  (a)  and  (b)  to  a 
nearly  isotropic  one  in  (c),  and  further  to  a  clearly  negative 
P  parameter  in  (d).  Since  a  two-photon  transition  is  invoked 
in  process  (2),  the  angular  distribution  was  analysed  by  the 
expression  [18-20] 

1(0)  oc  ^2^2 (cos  0)  +  /?4P4(cos0).  (3) 

Here  P2(cos  0)  and  P4(cos  0)  are  the  Legendre  polynomials  of 
degrees  2  and  4,  respectively,  fa  and  fa  are  the  corresponding 
anisotropy  coefficients  and  0  is  the  angle  between  the 
directions  of  the  laser  polarization  and  the  ejected  electron. 
The  best-fitted  distributions  (the  solid  lines)  along  with  the 
experimental  data  for  the  four  images  in  figure  3  are  displayed 
in  the  polar  representation  in  the  upper  panel  of  figure  5(a). 
The  lower  portion  of  figure  5  summarizes  the  best-fitted  fa 
and  fa  coefficients  over  the  spectral  profile  of  peak  6.  It  is 
quite  remarkable  that  the  photoelectron  angular  distributions, 
or  the  fa  and  fa  coefficients,  exhibit  such  dramatic  variations 
across  an  autoionization  resonance. 

Similar  measurements  and  analysis  were  performed  for 
peak  11,  and  the  final  results  are  presented  in  figure  5(b). 
Peaks  6  and  11  have  been  assigned  [13]  to  the  transitions 
from  2Pi/2  to  (3Pi)8p[2]5/2  and  to  (3Pi)8p[2]3/2,  respectively. 
(The  notation  of  the  superexcited  Rydberg  state  follows  the  Jcl 
coupling  scheme  [2 1  ] ,  in  which  the  total  angular  momentum  Jc 
of  the  I+  ion  core  is  coupled  to  the  orbital  angular  momentum 
/  of  the  Rydberg  electron  to  give  the  resultant  K.  The  spin, 
s,  of  the  Rydberg  electron  is  then  coupled  to  K  to  give 
the  total  angular  momentum  J.  The  state  is  thus  labelled  as 
(3P jc)nl[K]j .)  Therefore,  the  observed  autoionization  process 
(2)  corresponds  to  ejecting  the  8p  Rydberg  electron  to  form 
the  final  ion  state  I+(3P2).  Since  the  ion-core  state  changes 
from  3 Pi  to  3P2,  this  autoionization  process  must  invoke  the 
spin-orbit  couplings  of  the  3 Pi -core  Rydberg  states  to  the  3P2 
ion-continuum.  Physically,  it  can  arise  from  either  the  spin- 
flip  in  the  core  to  release  the  Rydberg  electron  or  the  exchange 
interactions  between  the  Rydberg  and  core  electrons  [22,  23]. 

Assuming  a  Coulomb  repulsive  form  for  the 
autoionization  interaction  matrix  and  approximating  it  as  a 
two-electron  problem,  a  propensity  rule  of  A /  =  0,  ±2  for 
autoionizing  electrons  has  previously  been  derived  [13].  For 
the  p-like  Rydberg  states,  such  as  the  present  peaks  6  and  11, 
the  dominant  outgoing  electron  waves  will  be  of  the  p  and 
f  types.  Hence,  the  observed  variations  of  the  photoelectron 
angular  distributions  across  the  autoionization  states  seen  in 
figure  5  arise  from  the  interferences  of  the  two  outgoing  waves 
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Figure  4.  The  upper  panel  presents  the  photoelectron  yield  spectrum  of  peak  6  and  the  fractional  branching  (•)  of  the  total  electrons  that  is 
from  autoionizing  the  I-atom.  The  lower  panel  shows  the  resultant  photoelectron  yield  spectrum  for  process  (2),  along  with  the 
photoionization  spectrum  of  the  I-atom  by  monitoring  I+.  For  comparison,  both  spectra  are  normalized  to  the  same  peak  heights. 


Figure  5.  Summary  of  the  autoionizing  lineshape  and  the  anisotropy  coefficients  of  photoelectron  angular  distribution  for  peak  6  in  (a)  and 
peak  1 1  in  (b).  The  points  are  the  experimental  data  and  the  lines  are  the  fits  (see  the  text  for  details).  Also  exemplified  in  the  upper  panels 
are  a  few  photoelectron  angular  distributions  in  the  polar  representation  at  different  energies  marked  by  arrows. 


whose  magnitudes  and  phases  are  continuously  changing 
as  the  photon  energy  changes.  Note  that  even  at  energies 
relatively  far  away  from  the  resonance  peak,  the  two-photon 
ejected  ‘background’  f-wave  could  conceivably  invoke  a 
shape  resonance  trapped  by  the  l  =  3  centrifugal  barrier — in 
contrast  to  the  Feshbach  type  associated  with  an  autoionizing 
resonance — so  that  the  rapid  variation  of  the  phase  shift  of  the 
direct  ionization  amplitude  might  need  to  be  considered.  It 
is  also  worth  noting  that  while  the  states  for  peaks  6  and  1 1 
invoke  the  same  Rydberg  orbital  and  the  same  K ,  the  degrees 
of  asymmetry  of  their  lineshapes  are  clearly  different,  so  as 
the  patterns  of  photoelectron  angular  distributions.  Inspection 


of  the  displayed  photoelectron  angular  distributions  suggests 
more  f-wave  characters  for  peak  11  than  peak  6.  Further 
theoretical  work  will  be  needed  to  understand  this  intriguing 
/-dependence  of  the  discrete-continuum  couplings. 

The  observed  asymmetric  lineshapes  (the  upper  panels, 
figure  5)  can  be  parametrized  using  the  Fano  profile  expression 
(equation  (1)).  The  resultant  fits  are  shown  as  the  solid  lines 
in  figure  5,  and  the  parameters  summarized  in  table  1.  Since 
this  is  a  multi -variable  fit  to  equation  (1),  the  choices  of  initial 
parameters  were  guided  by  the  extremes,  such  as  the  magnitude 
of  the  absorption  cross  section  far  from  the  peak  for  as 
well  the  peak  position/ width  and  the  degree  of  asymmetry 
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Table  1.  Fano  parameters  in  fitting  the  two-photon  absorption  cross  section,  equation  (1)  or  (6),  and  the  anisotropy  coefficients  of  the 
photoelectron  angular  distribution,  equation  (4)  or  (7).  See  the  text  for  details.  Note  that  the  ionization  potential  (IP)  of  the  iodine  atom  is 
84  295.1  cm-1  [21]  and  that  of  the  spin-orbit  excited  I(2Pi/2)  is  7603.2  cm-1  above  the  ground  state  I(2P3/2);  thus,  the  listed  autoionization 
states  range  from  about  97  cm-1  to  2666  cm-1  above  the  IP  for  peaks  1  and  14,  respectively. 


Parameters 

Peak  6 

Peak  1 1 

Peak  1/2 

Peak  13/14 

Er  (cm  l) 

85  803.98 

86296.85 

84  392.05/84  476.65 

86918.55/86  961.45 

T  (cm-1) 

24.7 

8.1 

4.95/4.4 

5.45/2.5 

o 

ob 

0.19 

0.12 

0.044 

0.14 

0.019 

0.0011 

0.0036/0.0023 

0.0051/0.0041 

q 

8.34 

-39.7 

-54/43 

-36/29.2 

A  (cm-1) 

-9.36 

+1.46 

+10.1/— 4.7 

+3. 44/ -1.03 

T  (cm-1) 

63.7 

31.5 

70.83/41.3 

36.51/12.4 

Pi 

ob' 

0.67 

0.58 

0.49 

0.41 

°a 

-1.11 

-0.29 

0.44/— 0.41 

0.32/— 0.39 

0.84 

0.62 

-0.5/0.28 

-0.43/0.21 

A  (cm-1) 

-9.36 

+1.46 

—3.8/— 73.5 

—2.4/— 33.1 

T  (cm-1) 

63.7 

31.5 

15/80 

15/50 

Pa 

Ob" 

0.11 

-0.65 

-0.2 

-0.41 

°a' 

0.16 

0.31 

0.27/— 0.001 

0.21/— 0.0025 

q" 

-0.46 

1.31 

1.11/11.6 

1.34/17.7 

for  q ,  as  illustrated  in  figure  1(b).  The  listed  parameters 
were  then  yielded  after  several  iterations  for  the  global  fit. 
The  uncertainties  for  Er  and  T  are  within  1  cm-1,  whereas 
the  three  parameters  cra ,  ob  and  q  are  highly  correlated  with  the 
estimated  uncertainty  of  ±2  for  the  q  value.  It  is  worth  noting 
that  from  the  fitted  widths  V,  the  effective  autoionization 
lifetimes  of  the  (3Pi)8p[2]5/2  (peak  6)  and  (3Pi)8p[2]3/2 
(peak  11)  Rydberg  states  can  be  estimated  to  be  0.22  ps  and 
0.65  ps,  respectively.  Compared  to  peak  6,  the  profile  of  peak 
1 1  appears  more  symmetric,  which  yields  a  larger  value  of  \q  | . 
Since  q  =  C\/ C2Ci2  and  from  table  1,  an  estimated  oa  (=C^) 
for  peak  1 1  is  an  order  of  magnitude  smaller  than  that  for  peak 
6;  a  larger  \q\  will  then  imply  a  significantly  smaller  value  of  C2 
and/or  Ci2.  However,  from  table  1  the  estimated  crb(=C%)  for 
the  two  peaks  are  comparable;  thus,  the  Rydberg-continuum 
coupling  (Ci2)  for  peak  1 1  must  be  weaker  than  that  for  peak 
6,  in  consistence  with  the  above  longer  lifetime. 

In  order  to  analyse  in  detail  the  measured  photoelectron 
angular  distribution,  we  followed  the  theoretical  prediction 
of  Grum-Grzhimailo  et  al  [12]  by  expressing  the  anisotropy 
coefficients  in  terms  of  a  Fano-like  form: 

ft  =of+ff0A(9A+e)2/(  1+g2).  (4) 

In  equation  (4)  the  parameters  a[’’ ,  a„  and  are  specific 
quantities  for  each  individual  observable  fo.  However,  the 
scaled  energy  s  =  (E  —  Er)/(T/2)  in  equation  (4)  is  fixed 
by  a  scaled  width  f  =  xV,  and  the  shifted  resonance  energy 
Er  =  Er  +  A  from  fitting  the  spectral  profile  (equation  (1)) 
with  *  =  [crb(cro  +  cfaq2)]l/2/o^  A  =  -(qTcra) /2cr0  and 
oo  =  cfa+ab  [12].  The  resultant  fits  to  the  energy  dependences 
of  the  anisotropy  coefficients  for  both  peaks  6  and  11  are 
displayed  as  the  solid  lines  in  the  respective  lower  portions 
of  figure  5,  and  the  parameters  are  also  listed  in  table  1.  As 
can  be  seen,  the  energy  dependences  of  both  /32  and  $4  are 
reasonably  reproduced,  despite  vastly  different  behaviours. 
More  remarkable  is  the  fact  that  the  two  parameters  f  and 
Er  used  in  fitting  both  fi2  and  @4  spectra  follow  exactly  the 


universal  scaling  relationships  given  above,  in  strong  support 
of  the  theoretical  prediction. 

An  alternative  approach  to  account  for  the  variation  in  p 
over  an  autoionizing  resonance  is  the  so-called  KS  formula 
[24,  25],  which  is  derived  for  the  one-photon  ionization 
processes.  In  that  framework 

P  =  —2(Xs2  +  Ys  +  Z)/(Ae2  +  Be  +  C),  (5) 

and  the  coefficients  A,  B  and  C  are  related  to  the 
Fano  parametrization  for  the  photoionization  cross  section, 
equation  (1),  as  A  =  oa  +  crb,  B  =  2 qoa  and  C  =  oaq2  +  ob. 
The  coefficients  A,  Y  and  Z  are  treated  as  free  parameters 
in  the  fit  to  the  data.  (Note  that  either  equation  (4)  or  (5) 
involves  three  fitting  parameters.)  While  we  found  that  the  KS 
formula  can  also  be  adopted  in  fitting  the  p2  spectra  for  peaks 
6  and  11,  its  theoretical  validity  to  fit  the  p4  spectrum  is  yet 
to  be  established.  Moreover,  the  use  of  equation  (5)  becomes 
ambiguous  in  the  case  of  overlapping  resonances  ( vide  infra). 
We  therefore  did  not  pursue  further  along  that  line. 

3.3.  Overlapping  resonances 

Figure  6  (the  upper  panels)  presents  two  pairs  of  autoionizing 
resonances,  peaks  1  and  2  in  (a)  and  peaks  13  and  14  in  (b). 
The  two  peaks  in  each  pair  appear  well  separated,  yet  not  far 
enough  to  allow  a  separated  fitting  to  each  individual  peak 
according  to  equation  (1).  To  simultaneously  fit  the  pair  of 
overlapping  resonances,  we  tentatively  used  a  Fano-like  form: 

a  =  ab  +oai(qx  +  ex)2  /  (l  +  s\)  +  cra2(q2  +  e2)2/ (l  +  e|) . 

(6) 

In  other  words,  we  assumed  that  the  direct  couplings  between 
the  two  resonance  states  are  negligibly  small,  yet  both 
resonances  are  coupled  to  a  common  continuum.  The  results 
are  shown  as  solid  lines  in  the  upper  portions  of  figure  6,  and 
the  fitting  parameters  are  listed  in  table  1.  From  the  fitted 
widths  T,  the  effective  lifetimes  of  those  autoionization  states 
range  from  about  2.1  ps  for  peak  14  to  0.97  ps  for  peak  13, 
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Figure  6.  Same  as  figure  5,  but  for  two  pairs  of  overlapping  resonance  states:  (a)  for  the  pair  of  peaks  1  and  2,  and  (b)  for  peaks  13  and  14. 
Note  that  the  width  and  shift  in  fitting  the  fa  spectra  follow  the  scaling  relationships,  but  they  are  treated  as  free  parameters  for  the  fitted  fa 
spectra. 


i.e.  somewhat  longer  than  the  above  isolated  resonances.  More 
rigorous  treatment  will  have  to  be  rendered  to  theoretical  works 
on  overlapping  resonances  [26,  27]. 

The  energy  dependences  of  the  corresponding  anisotropy 
coefficients  are  shown  in  the  lower  panels.  An  analogous  form 
is  assumed  in  fitting  the  /3  spectra,  i.e. 

Pi  =  Ob  +  aa\  +  gl)2/ (l  +  gl)  +  °a2  (<l2  +  g2)2/ (l  +  el)- 

(7) 

Note  that  the  universal  scaling  relationships  proposed  by 
Grum-Grzhimailo  are  formulated  for  an  isolated  resonance; 
it  is  not  clear  a  prior  if  they  will  also  hold  for  overlapping 
resonances.  Nevertheless,  using  the  universal  scaling 
relationships,  we  found  that  equation  (7)  can  reproduce  the 
energy  dependences  of  the  fa  parameters  for  both  pairs 
(shown  as  the  solid  lines  in  the  lower  portions  of  figure  6), 
but  it  failed  to  account  for  the  fa  spectra — in  particular  for 
the  peculiar,  central  peaking  behaviour  around  the  resonance 
energy  of  peaks  1  and  13,  respectively.  However,  if  the  shifts 
and  the  widths  for  fa  are  treated  as  free  parameters,  the 
energy  dependences  of  fa  can  be  reasonably  reproduced  by 
equation  (7),  as  shown  by  the  solid  lines  in  figure  6. 

Two  plausible  reasons  can  be  offered  for  the  failure  of 
the  universal  scaling  relationship  in  fitting  the  fa  parameters. 
Firstly,  the  expression  of  equation  (6)  or  equation  (7)  assumes 
a  case  of  two  non-interacting  resonance  states  coupled  to 
one  continuum.  Hence,  its  form  may  not  be  applicable  if 
the  coupling  between  the  two  overlapping  resonances  cannot 
be  neglected.  Secondly,  the  autoionization  states  for  peaks 
1  and  2  are  (3Pi)4f[2]3/2,5/2  and  (3Pi)4f[3]5/2,  respectively 
[13,  21].  In  other  words,  accounting  for  the  two  fine-structure 
states  hidden  underneath  peak  1,  three  resonance  states  are 
actually  involved — rather  than  two  resonances  as  implicated  in 
equations  (6)  and  (7).  Two  of  them  are,  however,  not  resolvable 
in  the  photoelectron  spectra  (i.e.  shown  as  peak  1)  nor  can  they 
manifest  themselves  as  any  noticeable  feature  in  the  energy 


dependences  of  the  fa  parameters,  which  suggests  that  the 
direct  coupling  between  the  two  heavily  overlapped  resonances 
may  need  to  be  considered,  rendering  the  mere  addition  of 
a  third  resonance  term  in  equations  (6)  and  (7)  inadequate. 
Further  theoretical  works  will  be  required.  Nevertheless, 
the  observation  of  an  intriguing  local  structure  near  peak  1 
in  the  energy  dependence  of  fa  seems  to  suggest  this  as  a 
more  likely  explanation.  Based  on  the  previously  derived 
propensity  rule  A /  =  0,  ±2  for  autoionization  [13],  the  f- 
like  Rydberg  state  will  preferentially  eject  p-,  f-  and  h- waves. 
Thus,  three  outgoing  waves  are  involved  in  the  vicinity  of  each 
autoionization  resonance  state,  and  nine  waves  in  total  need 
to  be  considered  over  the  energy  range  of  figure  6.  More 
perplexing  interference  patterns  are  anticipated.  Inspection 
of  the  exemplified  angular  distributions  in  the  upper  panel  of 
figure  6(a)  indicates  that  it  is  indeed  the  case.  Similarly,  the 
autoionization  states  for  peaks  13  and  14  are  (3Pi)5f[2]3/2,5/2 
and  (3Pi)5f[3]5/2,  be.  the  same  respective  Rydberg  series  as 
peaks  1  and  2  [  1 3] .  It  is  therefore  not  too  surprising  to  note  that 
the  general  patterns  of  fa  and  fa  spectra,  or  the  overall  shapes 
of  the  positive  fa  and  negative  fa  values,  are  remarkably 
similar  for  the  pair  of  peaks  1  and  2  and  peaks  13  and  14. 

4.  Summary 

An  in-depth  study  of  autoionizing  dynamics  of  the  iodine 
atom  by  two-photon  excitation  was  carried  out  using  the 
photoelectron  imaging  technique.  Particular  emphasis  is 
placed  on  the  energy  dependences  of  the  ejected  photoelectron 
angular  distributions.  A  dramatic  variation  was  found  in 
the  vicinity  of  autoionization  resonance  states.  Exploiting 
the  scaling  relationship  predicted  theoretically,  the  energy 
evolutions  of  the  observed  anisotropy  coefficients  can  be 
accounted  for  in  the  case  of  isolated  resonances.  However, 
in  the  case  of  overlapping  resonances,  while  the  scaling 
relationship  remains  held  for  the  fa  parameters,  it  fails  to 
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reproduce  the  fa  behaviours.  Possible  explanations  are 
suggested,  and  further  theoretical  investigations  are  warranted. 
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The  vibrational  branching  ratio  [CH3(v  =  0)  +  HCl(v'  =  1)]/ 
[CH3(v  =  0)  +  HCl(v'  =  0)]  of  two  correlated  product  pairs 
from  the  title  reaction  shows  a  dramatic  /^-dependence,  in  sharp 
contrast  to  the  previously  observed  behavior  in  Cl  +  CHD3- 
(vx  =  1),  while  the  vibrational  enhancement  factors  in  reactivity 
of  the  two  reactions  are  remarkably  similar. 

It  has  been  well  documented,  both  experimentally1  11  and 
theoretically,12  18  that  excitation  of  different  vibrational  mo¬ 
tions  of  a  polyatomic  reactant  can  exert  profound  effects  on 
chemical  reactivity.  What  emerged  from  those  mode-selective 
studies  is  the  intuitively  appealing  picture  that,  in  a  direct 
bimolecular  reaction,  the  vibrational  state  containing  large 
components  of  motion  along  a  particular  reaction  coordinate 
preferentially  drives  the  system  over  the  transition  state  along 
that  coordinate,  leading  to  mode-dependent  reactivity  and 
product  branching  ratios.  Among  the  systems  that  have  been 
studied,  the  reaction  of  chlorine  atom  with  methane  and  its 
isotopic  variants  stands  as  a  benchmark.5  16  In  the  context  of 
mode-specific  behaviors,  the  reaction  with  CH4  is  of  interest  in 
that  CH4  has  two  nearly  degenerate  stretching  modes,  the 
symmetric  (vi)  and  the  antisymmetric  (v3),  and  it  is  not  totally 
clear,  a  priori ,  which  mode  will  be  more  effective  in  promoting 
the  reaction  rate.  At  the  outset,  we  should  point  out  that  the 
vibrational  eigenstates  are  named  by  their  dominant  zero- 
order  state  characters.  The  nominal  v\  and  v3  fundamental 
stretching  modes  of  CH4  contain  small  admixture  of  bending 
mode  character  (~4-8%  of  contributions)  in  higher  order 
polyads. 19,20 

In  a  pioneering  work,5  Zare  and  coworkers  examined  the 
effect  of  antisymmetric  stretch  excitation  of  CH4(v3  =  1)  on 
the  reaction  cross  section  and  scattering  dynamics  and  ob¬ 
tained  a  vibrational  enhancement  of  30  ±  15  at  a  mean 
collision  energy  Ec  =  15.5  kJ  mol-1,  which  was  attributed 
to  the  opening  of  the  cone  of  acceptance  of  the  reaction  upon 
vibrational  excitation.21  A  very  recent  re-investigation  refined 
the  enhancement  factor  to  13.4  ±  2. 4. 22  A  later  study  from  the 
same  laboratory  compared  the  dynamical  effects  of  symmetric 
(vi)  and  antisymmetric  (v3)  stretch  excitations  of  CH4  in  Cl  + 
CH4.6  They  found  that  both  stretching  excitations  yield  vir¬ 
tually  identical  product  internal  state  and  angular  distribu¬ 
tions.  Based  on  these  observations,  they  concluded  that, 
whereas  the  initial  preparation  of  CH4  in  different  stretching 
modes  might  alter  the  total  reactivity  (vide  infra),  the  dynamics 
leading  to  product  formation  follow  a  common  pathway, 
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which  could  be  promoted  by  vibrational  mixing  during  the 
collision  event,  scrambling  the  initial  mode-selection. 

In  that  study,  the  infrared  (IR)-inactive  v\  =  0  -►  1 
excitation  of  CH4  was  achieved  by  stimulated  Raman  pump¬ 
ing.  Comparison  of  its  reactivity  enhancement  to  the  v3- 
excitation  could  not  be  made  because  the  fraction  of  methane 
molecules  being  pumped  was  not  quantified.  The  difficulty  of 
comparing  the  relative  reactivities  of  the  two  fundamental 
stretching  modes  was  alleviated  by  Crim  and  coworkers  in 
studying  the  two  stretch-bend  combination  vibrations  of  CH4, 
v\  +  v4  and  v3  +  v4  (v4  is  the  umbrella  bending  vibration),  for 
which  both  modes  can  be  activated  by  direct  IR  excitations.7 
They  observed  a  noticeable  difference  in  reactivities:  cr(vi  +  v4) 
=  1.9  x  o-(v3  +  v4)  =  19  x  o(v  =  0),  arguing  different 
reaction  paths  for  the  two  types  of  vibrational  excitations. 
They  further  proposed  a  vibrationally  adiabatic  model  to 
rationalize  the  differential  reactivities  of  the  two  modes  of 
excitation.7,18  In  this  model,  the  approach  of  the  Cl-atom 
causes  the  vibrational  energy  of  the  CH4  symmetric  and 
antisymmetric  stretches  to  become  localized  in  the  proximal 
(reactive)  and  distal  (unreactive)  C-H  bonds,  respectively. 
Thanks  to  this  collision-induced  intramolecular  vibration  en¬ 
ergy  redistribution  (IVR),  the  symmetric  stretch  has  more 
energy  along  the  reaction  coordinate,  leading  to  a  higher 
reactivity  than  the  antisymmetric  stretch  excitation.  In  other 
words,  the  reaction  remains  largely  vibrationally  adiabatic 
without  significant  state  mixing  of  two  local  modes  during 
reactive  collisions.  Since  there  may  be  a  synergic  effect  when  a 
combination  mode  is  excited,  the  subtle  differences  between 
the  above  two  viewpoints  (i.e.,  same  reaction  path  or  not  and 
the  nature  of  mode  mixings)  remain  to  be  reconciled. 

Another  motivation  for  initiating  the  present  study  is  from 
one  of  our  earlier  works  on  the  ground-state  reaction.23 
Examining  the  /^-evolution  of  two  product  pair-correlated 
angular  distributions  over  a  wide  energy  range  from  about  8.4 
kJ  mol-1  to  84  kJ  mol-1  revealed  vastly  distinct  patterns. 
While  the  result  for  the  CH3(v  =  0)  +  HCl(v'  =  0)  pair  shows 
a  pronounced  ridge- structure  suggestive  of  peripheral  dy¬ 
namics  for  a  reaction  governed  by  large  impact-parameter 
collisions,24,25  that  for  the  CH3(v  =  0)  +  HCl(v'  =  1)  pair 
exhibits  a  pattern  characteristic  of  resonance  reaction  mechan¬ 
ism.26  Based  on  previous  ab  initio  calculations,12,13  the  corre¬ 
lation  between  the  reactant  and  product  vibrational  modes 
was  then  constructed.11,23  A  dynamic  potential  well  correlat¬ 
ing  the  reactant  Cl  +  CH4(v!  =  1)  to  the  CH3(v  =  0)  + 
HCl(v'  =  1)  product  pair  was  found,  in  support  of  the 
proposed  reactive  resonance  mechanism.  The  reaction  of 
Cl  +  CH4(v  =  0)  is,  however,  largely  vibrationally  adiabatic 
even  at  Ec  ~  84  kJ  mol-1;  namely,  more  than  98%  of  reactive 
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fluxes  proceed  through  the  vibrational  ground-state  reaction 
path  leading  to  the  CH3(v  =  0)  +  HC1(V  =  0)  product  pair.23 
Thus,  merely  2%  of  reactive  trajectories  are  mediated  by 
reactive  resonance.  To  shed  more  light  on  the  proposed 
reactive  resonance,  it  is  then  desirable  to  initiate  the  reaction 
with  CH4(v!  =  1)  for  better  access  of  the  dynamic  well. 
However,  an  efficient  Raman  pumping  for  v3  =  1  excitation 
is  more  demanding.  If  the  aforementioned  collision-induced 
IVR7,18  or  the  mode-mixing  mechanism6  is  at  work,  an 
investigation  with  direct  IR-pumped  CH4(v3  =  1)  could  then 
be  a  viable  alternative. 

Reported  here  is  the  first  product  pair-correlation  measure¬ 
ment  of  the  Cl  +  CH4(v3  =  1)  -►  CH3(v  =  0)  +  HCl(v'  = 
0,1)  reaction.  The  ground-state  reaction  of  Cl  +  CH4  is 
endothermic  by  5.1  kJ  mol-1,  with  a  reaction  barrier  about 
10.5  kJ  mol-1,25  whereas  one  quantum  of  excitation  in  the  v3 
and  v3  modes  provides  CH4  molecule  with  energy  equivalent 
to  34.9  kJ  mol-1  and  36.1  kJ  mol-1,  respectively.  To  elucidate 
the  mode-specific  behavior  of  this  reaction,  comparisons  were 
made  to  both  the  ground-state  reaction23  and  the  analogous 
reaction  of  Cl-atom  with  C-H  stretch-excited  CHD3(vi  = 
l).11  The  latter  comparison  serves  as  a  conceptually  interesting 
contrast  in  that  the  initially  prepared  CH4(v3  =  1)  is  a 
delocalized  normal-mode  vibration,  whereas  for  CHD3(vi  = 
1)  the  vibrational  energy  is  directly  deposited  into  the  localized 
C-H  bond  that  is  to  be  broken.  These  close  scrutinies  reveal 
rather  different  behaviors,  which  in  turn  lead  to  a  new 
perspective  about  the  effects  of  the  two  different 
stretching  modes. 

The  experimental  details  are  the  same  as  the  previous  study 
on  Cl  +  CHD3(vi  =  l),11  only  a  brief  description  is  given 
here.  A  discharge-generated,  pulsed  Cl-beam  (5%  Cl2  in  He  at 
6  atm  total  pressure)  was  double-skimmed  before  crossing 
with  a  CH4  beam  pulsed  from  another  valve.  A  tunable  IR 
optical  parametric  oscillator/amplifier  (OPO/A),  operated  in 
near  single-mode,  prepared  the  CH4  reactant  prior  to  the 
collision  center  with  one-quantum  excitation  of  antisymmetric 
stretch  via  the  v3  =  0  ->  1,  R(l)  transition  at  3038.50  cm-1.27 
The  absolute  frequency  of  the  OPO/A  output  was  calibrated 
with  the  aid  of  a  photoacoustic  cell  filled  with  ~70  torr  CH4. 
After  the  collision,  the  dominant  CH3(v  =  0)  product  was 
probed  by  a  (2+1)  resonance-enhanced  multiphoton  ioniza¬ 
tion  (REMPI)  process  at  333.57  nm  via  the  0°-band  of  the 
X2A2"  — ►  3p2A2//  Rydberg  transition,28  and  a  time-sliced  ion 
velocity  imaging  technique29  mapped  the  state-correlation  of 
coincidently  formed  HC1  coproducts.11,30,31  The  initial  colli¬ 
sion  energy  Ec  was  controlled  by  varying  the  intersection 
angles  of  the  two  molecular  beams.  To  interrogate  the  effects 
of  reactant  vibrational  excitation,  at  each  Ec  the  product 
images  were  acquired  with  IR  laser  on  and  off  alternatively 
to  minimize  the  long-term  drifts  and  other  possible  systematic 
errors. 

Fig.  1(a)  and  (b)  exemplify  two  raw  images,  IR-on  and  IR- 
off,  at  Ec  =  24.3  kJ  mol-1.  The  IR-off  image  is  virtually 
identical  to  that  acquired  previously:23,25  the  sideway-domi¬ 
nant  ringlike  feature  is  the  CH3(v  =  0)  +  HCl(v'  =  0)  product 
pair  with  small  contribution,  along  the  outer  rim,  of  the  same 
pair  from  reaction  with  bend-excited  CH4.  (Supersonic  expan¬ 
sion  yields  a  translationally  and  rotationally  cold  (~5-10  K) 


Fig.  1  Time-sliced  raw  images  of  the  probed  CH 3(z/  =  0)  products 
from  the  Cl  +  CH4  reaction  at  Ec  =  24.3  kJ  mol-1  for  IR-on  (a)  and 
IR-off  (b).  The  0°-angle  is  the  initial  direction  of  the  CH4  beam  in  the 
center-of-mass  (cm)  frame.  The  ringlike  features  of  the  product  images 
are  assigned  as  indicated.  Some  backgrounds  such  as  the  bright  spot  in 
the  backward  direction  can  readily  be  identified  and  were  discounted 
in  the  data  analysis,  (c)  The  angular  distributions  of  the  inner  ring 
structures  from  IR-on  (a)  and  IR-off  (b),  respectively,  (d)  The  derived 
angular  distributions  of  the  (0o,  1)#  and  (0o,  0)#  product  pairs  from 
the  Cl  +  CH4(j/3  =  1)  reaction. 


methane  beam,  but  with  thermally  populated  vibration  mani¬ 
folds.10)  When  the  IR  laser  is  turned  on,  two  new  features  are 
vividly  displayed.  On  energetic  grounds,  the  outer  ring  is 
ascribed  to  the  Cl  +  CH4(v3  =  1)  ->  CH3(v  =  0)  + 
HCfiV  =  0)  reaction  (denoted  as  the  (00,  0)^  product  pair 
with  the  supercript  #  indicating  the  stretch-excited  reactant) 
and  the  intense,  sharp  forward  peak  in  the  inner  ring  is  from 
the  CH3(v  =  0)  +  HCl(v'  =  1)  or  the  (00,  \)^  pair,  which  is 
nearly  degenerate  to  the  (00,  0)  pair  from  the  ground- state 
reaction  (IR-off  image). 

To  disentangle  the  overlapped  contributions  to  the  inner 
ring  of  the  IR-on  image,  we  followed  the  approach  recently 
developed  in  the  studies  of  the  Cl  +  CHD3(vi  =  1)  reac¬ 
tion.11,32  Using  the  threshold  method,32  we  obtained  about 
19%  of  CH4  reactants  being  stretch-excited.  By  scaling  down 
the  IR-off  angular  distribution  by  0.19  and  subtracting  it  from 
the  IR-on  data  (Fig.  1(c)),  the  genuine  distribution  for  (00,  1)^ 
was  then  recovered  from  the  overlapped  inner  ring.  The 
resultant  pair-correlated  angular  distributions  are  presented 
in  Fig.  1(d),  which  show  vastly  different  appearances  from  the 
(00,  0)^  pair,  strongly  suggesting  different  reaction  pathways. 
A  broad  back/sideways  angular  distribution  signifies  a  direct 
rebound  mechanism25  in  forming  the  (00,  0)^  product  pair  at 
this  collisional  energy,  whereas  a  sharply  forward  peaking 
distribution  for  the  (00,  1)^  pair  generally  suggests  time-delay 
reaction  mechanisms33  and  most  likely,  in  the  present  case,  a 
short-lived  complex  pathway  as  the  (00,  1)  pair  from  ground- 
state  reaction.23  Similar  measurements  were  performed  at 
different  Ec ,  and  four  of  them  are  presented  in  Fig.  2.  The 
overall  behaviors  are  quite  similar,  except  the  subtle  variation 
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Fig.  2  As  Fig.  1(d),  except  at  four  different  collision  energies,  Ec  in 
the  unit  of  kJ  mol-1. 


in  the  relative  intensities  of  the  sharp  forward  peak  to  the  more 
broadly  scattered  component  of  the  (0o,  \)^  pair.  It  is  worth 
noting  that  these  pair-correlated  angular  distributions  bear  a 
remarkable  resemblance  to  the  corresponding  ones  in  the  Cl  + 
CHD3(V!  =  1)  ->  CD3(v  =  0)  +  HC1(V  =  0,1)  reaction.1132 

Integrating  each  distribution  shown  in  Fig.  1  and  2  over  all 
angles,  weighted  by  the  sin 9  term  for  the  solid-angle  factor 
(where  6  is  the  product- scattering  angle),  and  accounting  for 
only  ~  19%  of  the  CH4  reactants  being  pumped,  we  deduced 
the  normalized,  pair-correlated  integral  cross  section  at  each 
Ec.  The  reactive  excitation  function  (J0(EC)  i.e.,  the  dependence 
of  the  integral  cross  section  on  Ec ,  of  the  ground-state 
reaction  obtained  from  the  IR-off  images  agreed  with  the 
previous  report.25  By  measuring  the  relative  cross  section  for 
CH4(v3  =  1)  to  that  for  CH4(v  =  0)  at  each  Ec,  the  reactive 
excitation  function  o^(Ec )  of  the  Cl  +  CH4(v3  =1)  reaction 
was  then  determined  and  presented  in  Fig.  3(a).  As  is  seen, 
while  the  ground-state  reaction  displays  a  clear  threshold 
around  10.5  kJ  mol-1,25  the  stretch-excited  reactivity  remains 
significant  over  the  ifc-range  of  this  study.  Plotted  in  Fig.  3(b) 
on  logarithmic  scale  is  the  vibrational  enhancement 
factor  (cr^/do)  which,  due  to  the  threshold  effect,  increases 
sharply  from  less  than  3  to  about  60  with  the  decrease  in  Ec.  At 
Ec  =  15.5  kJ  mol-1,  a  vibrational  enhancement  factor  of  ~  10 
is  in  reasonable  agreement  with  the  recent  result  of  13.4  by 
Zare  and  coworkers.22  It  is  also  worth  noting  that  the 
independence  of  the  vibrational  enhancement  factor  resem¬ 
bles  that  found  previously  for  the  Cl  +  CF[D3(v1  =  1) 
reaction,11  except  that  the  whole  curve  in  the  present  case 
appears  shifted  downward  by  a  factor  of  1.4  ±  0.2.  Clearly,  the 
interpretation  of  the  observed  vibrational  enhancement  in 
reactivity  goes  beyond  the  simple  collision-induced  IVR  me¬ 
chanism,7,18  in  which  the  initially  delocalized  v3-mode  excita¬ 
tion  will  be  transformed  into  a  predominantly  distal  or 
unreactive  local  mode.  It,  thus,  can’t  account  for  such  a 
significant  enhancement,  unless  subsequent  nonadiabatic  tran¬ 
sitions  occur  between  the  two  local  modes. 

The  comparison  of  the  two  isotopic  reactions  is 
more  striking  in  terms  of  pair-correlated  vibration  branching 
erf  /cr0,  where  =  Gq  +  erf  and  erf  denotes  the  reactive 


Fig.  3  (a)  Normalized  reactive  excitation  functions  of  the  Cl  + 
CH4(i/)  ->  CH3(^  =  0)  +  HCl(^)  reactions,  with  o #  and  cr0  denoting 
the  reactive  cross  sections  from  the  antisymmetric-stretch  excited  (z/3 
=  1)  and  ground-state  ( v  =  0)  reactants,  respectively,  (b)  The 
vibrational  enhancement  factors  at  fixed  Ec.  Also,  shown  by  “x”, 
for  comparison  is  the  recent  measurement  by  Zare  and  coworkers 
(ref.  22)  at  Ec  =  15.5  kJ  mol-1,  (c)  The  vibrational  branching  fractions 
of  the  coincidently  formed  HC1  {y'  =  1)  products  from  the 

Cl  +  CH4(^3  =  1)  reaction.  Solid  lines  are  visual  guides. 

cross  section  in  forming  the  CH3(v  =  0)  +  HCl(v'  =  i) 
product  pair  from  the  stretch-excited  reactant.  While  the 
vibrational  branching  fraction  for  the  Cl  +  CHD3(vi  =  1) 
— ►  CD3(v  =  0)  +  HCl(v')  reaction  rises  abruptly  at  energetic 
threshold  and  then  remains  roughly  constant  around 
0.44, 11  the  present  reaction  (Fig.  3(c))  displays  a  clear  down¬ 
ward  trend  from  about  0.65  at  Fc  —  10  kJ  mol1  to  0.4 
at  Ec  ~  28  kJ  mol-1. 

Based  on  the  vibrational  mode  correlation  between  the 
reactant  and  product  pairs,  we  previously  proposed  a  con¬ 
ceptual  framework  for  understanding  mode-specific  reactiv¬ 
ity.11  In  this  framework,  the  reaction  of  Cl  +  CHD3(vi  =  1) 
adiabatically  correlates  to  the  CD3(v  =  0)  +  HCl(v'  =  1) 
product  pair.  Hence,  the  formation  of  the  ground-state  pro¬ 
duct  pair  CD3(v  =  0)  +  HCl(v'  =  0),  which  correlates  to  the 
ground-state  reactants,  from  the  CHD3(v3  =  1)  reaction  must 
involve  vibrationally  nonadiabatic  transitions,  in  line  with 
theoretical  predictions  of  strong  curvature  couplings  between 
the  C-H  stretching  (v3  =  1)  mode  and  the  reaction  coordi¬ 
nate.12  14  In  other  words,  the  curvature  coupling  of  the  C-H 
stretching  mode  induces  a  bifurcation  of  the  reactive  fluxes 
into  an  adiabatic  and  a  nonadiabatic  pathways  leading  to  the 
observed  (00,  1)^  and  (00,  0)^  product  pairs,  respectively.11 
Nearly  constant  vibration  branching  fraction  erf  jo^  =  0.44 
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over  the  energy  range  of  10-50  kJ  mol  1  will  then  suggest  that 
the  curvature  coupling  is  not  a  strong  function  of  Ec. 

By  contrast,  the  present  reaction  of  Cl-atom  with  a  triply 
degenerate,  antisymmetrically  stretch-excited  CH4(v3  =  1) 
correlates  adiabatically  to  the  stretch-excited  CH3(v!  =  1 
and  v3  =  1)  and  HC1(V  =  0)  product  pair.23  Ab  initio 
calculations  showed  that  the  generalized  vibration  frequencies 
of  that  correlation  are  nearly  invariant  along  the  reaction 
coordinate.12,13  It  implies  that  the  reactant  v3-mode  is  a 
spectator;  thus,  no  vibrational  enhancement  is  anticipated 
for  that  adiabatic-correlated  product  channel.  Here,  the 
probed  product  pair  CH3(v  =  0)  +  HC1(V  =  1)  correlates 
instead  to  the  Cl  +  CH4(vi  =  1)  reactant  pair,  for  which  an 
adiabatic  potential  well  is  present  in  the  transition  state 
region.23  As  the  Cl-atom  approaches  the  IR-excited  CH4(v3 
=  1),  the  intermolecular  interactions  will  perturb  the  vibra¬ 
tional  motions  of  an  isolated  molecule.  As  such,  the  laser- 
prepared  eigenstate  evolves  into  a  superposition  state  in  the 
basis  of  vibrational  eigenstates  for  the  C1-CH4  colliding  pair, 
and  mode-mixing  ensues.  In  other  words,  the  formations  of 
(00,  1)^  and  (00,  0)^  product  pairs  from  the  v3-excited  reactant 
are  conceivably  initiated  by  the  vi/v3  mode-mixing  in  the 
entrance  valley  of  the  Cl  +  CH4  reaction,  as  previously 
suggested.6,7  This  collision-induced  mode  mixing  in  the  en¬ 
trance  valley  is  expected  to  be  facile  for  both  the  proximity  of 
the  vibrationally  adiabatic  surfaces  of  the  two  stretches  and 
long  interaction  time  of  heavy  reactants’  slow  motions.  Since 
the  excited  vr -adiabatic  pathway  has  a  lower  activation  bar- 
rier  and  been  predicted  theoretically  more  reactive,  a 
substantial  mixing  of  the  vr -character  into  the  initially  pre¬ 
pared  v3-excited  reactant  could  then  result  in  the  rate  promo¬ 
tion  of  the  Cl  +  CH4(v3  =1)  reaction  as  seen  in  Fig.  3(b). 

However,  the  observation  of  a  relatively  high  branching 
fraction  in  forming  the  (00,  1 pair  near  the  energetic  thres¬ 
hold  followed  by  a  rapid  decline  with  the  increase  in  Ec 
(Fig.  3(c))  is  perplexing  and  in  sharp  contrast  to  the  behaviors 
observed  in  the  Cl  +  CHD3(v!  =  1)  reaction.11  This  un¬ 
expected  result  suggests  that  the  reactivity  of  CH4(v3  =  1) 
leading  to  different  product  pairs  is  probably  more  than  just 
the  vi/v3  mode  mixing;  otherwise,  an  isc-independent  vibration 
branching  as  that  from  the  C-H  stretch-excited  CHD3  reac¬ 
tion11  would  have  been  observed.  Additional  factors  must  be 
involved  in  differentiating  the  vibrational  branching  behaviors 
of  the  two  stretch-excited  isotope  reactions. 

It  is  also  instructive  to  compare  the  HC1  vibrational  branch¬ 
ing  fractions  from  the  stretch-excited  and  ground-state  reac¬ 
tants  of  the  two  isotopic  reactions.  Despite  the  different 
dependences  on  Ec  as  just  discussed,  the  correlated  vibration 
branchings  from  both  stretch-excited  reactants  are  more  than 
one  order-of-magnitude  larger  than  the  respective,  transla- 
tionally  hot  ground-state  reactions  with  equivalent  amount  of 
total  energy.11,23  For  the  ground-state  reaction,  based  on  the 
observed  characteristic  pattern  on  the  dcr/d(cosd)  —  0  —  Ec 
plot  of  the  (00,  1)  product  pair,  a  resonance-complex  reaction 
mechanism  was  proposed  for  its  formation.23  And  the  path¬ 
way  leading  to  the  resonant  complex  formation  is  induced  by  a 
translation-to- vibration  (T  ->  V)  energy  transfer  process  in  the 
entrance  valley.23  If  the  behavior  of  CH4(v3  =  1)  reactivity  is 
indeed  mediated  in  part  (i.e.,  although  it  can’t  be  the  full  story 


as  aforementioned)  by  the  vi/v3  mode-mixing,  then  the  large 
disparity  in  vibrational  branching  fractions  between  the  two 
methane  states  (the  antisymmetric  stretch-excited  and  ground 
states)  may  reflect  the  relative  inefficiency  of  the  T  ->  V  energy 
transfer  compared  to  the  collision-induced  mode  mixing  in  the 
entrance  valley. 

In  summary,  product  pair-correlation  was  measured,  using 
a  time-sliced  velocity-map  imaging  technique,  for  the  Cl  + 
CH4(v3  =  1)  -►  CH3(v  =  0)  +  HCl(v'  =  0,1)  reaction  under 
crossed-beam  conditions.  Over  the  energy  range  of  this  study, 
significant  vibrational  enhancements  in  reactivity  were  ob¬ 
served.  Its  energy  dependence  is  nearly  the  same  as  that  for 
the  isotopically  analogous  reaction  of  Cl  +  CHD3(vi  =  1)  — ► 
CD3(v  =  0)  +  HCl(v'  =  0,1)  except  an  overall  downward 
shift  by  a  factor  of  1.4.  The  correlated  HC1  vibrational 
branching  fraction  erf /er^,  however,  exhibits  an  intriguing 
F'c-dependency,  in  sharp  contrast  to  the  E^-independent 
branching  for  the  Cl  +  CHD3(vi  =  1)  reaction.  The  origin 
of  this  difference  is  not  fully  understood  at  present.  Further 
study  to  unveil  the  nature  of  state  mixings  is  in  progress. 
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By  using  a  time-sliced  (ion)  velocity-imaging  technique,  the  title  reaction  was  investigated  to 
interrogate  the  pair-correlated  dynamics  of  HC1  (DC1)  in  concomitance  with  the  CHD2  or  CHD2 
ground-state  product  under  the  crossed-beam  conditions.  Product  pair-correlated  excitation 
functions,  vibrational  branching  ratios,  and  angular  distributions  were  obtained  over  a  wide  range 
of  collision  energies,  from  2  to  22  kcal/mol.  Two  distinct  reaction  mechanisms  were  uncovered. 
Both  the  dominant  ground-state  reaction  of  Cl+CH2D2(u  =  0)  -^HCl(u  =  0)/DCl(i;  =  0) 
+  CHD2(00)/CH2D(00)  and  the  reaction  forming  the  same  product  pairs  from  bend-excited  CH2D2 
reactants  proceed  through  direct  mechanisms,  which  shift  from  a  rebound  dynamics  near  thresholds 
to  a  peripheral  dynamics  at  higher  Ec.  Reactivity  of  forming  the  HCl(u  =  l)/DCl(u  =  1 ,2) 
+  CHD2(00)/CH2D(00)  product  pairs  from  ground-state  reactants  is  small,  yet  shows  dynamics 
pattern  characteristics  of  a  resonant  reaction  pathway.  Significant  spin-orbit  reactivity  of  Cl* (2P1/2) 
was  discovered,  and  its  mechanism  appears  to  be  also  mediated  by  the  resonant  reaction  pathway. 
Comparing  the  dynamical  attributes  of  the  H-  and  D-atom  transfer  channels,  remarkable  isotope 
effects  are  found — in  line  with  the  previous  findings  in  the  C1  +  CH4/CD4  and  C1  +  CHD3 
reactions.  ©  2008  American  Institute  of  Physics.  [DOI:  10.1063/1.2993264] 


I.  INTRODUCTION 

The  hydrogen  abstraction  reaction  of  the  chlorine  atom 
with  methane  plays  a  crucial  role  in  the  ozone  production/ 
depletion  cycle  in  the  stratosphere1  and  has  been  the  subject 
of  intense  investigations,  both  experimentally2  and 
theoretically,3  over  the  past  decades.  More  recently,  this  re¬ 
action  has  been  under  a  more  detailed  dynamical  scrutiny 
and  is  becoming  a  paradigm  for  fundamental  understanding 
of  polyatomic  reaction  dynamics. 4-6  In  particular,  ample  evi¬ 
dence  has  been  gathered,  demonstrating  the  dominance  of 
the  peripheral  reaction  mechanism.  A  peripheral  mecha¬ 
nism  occurs  for  reactions  proceeding  preferentially  with 
large  impact-parameter  ( b )  collisions,  giving  rise  to  a  skewed 
opacity  function  P(b)  strongly  in  favor  of  higher  b  values.7 
In  order  for  a  large  impact-parameter  trajectory  to  be  reac¬ 
tive,  the  attacking  Cl  atom  needs  not  only  to  move  under  a 
central  potential  but  also  to  experience  an  attraction  to  the 
abstracted  peripheral  hydrogen  atom. 

One  of  the  experimental  manifestations  of  peripheral  dy¬ 
namics  appears  as  the  unusual  shape  of  the  product  angular 
distribution.7-11  Although  the  distribution  is  mainly  back¬ 
ward  scattered  near  the  reaction  threshold  for  a  limited  range 
of  small  impact-parameter  collisions  to  reaction,  as  the  col¬ 
lision  energy  Ec  increases,  it  shifts  progressively  toward 
sideways/forward  with  a  characteristic  sharp  cutoff  in  angles 
against  the  formation  of  small  scattering-angle  products.  In 
the  primitive  line-of-center  model,  this  sharp  cutoff  origi¬ 
nates  from  the  requirement  that  only  those  collisions  with 
sufficient  line-of-center  energy  to  surmount  the  centrifugal- 
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shifted  barrier  can  react,8  and  the  characteristic  angular  dis¬ 
tribution  with  predominant  sideways  peaking  arises  from  the 
one-to-one  correspondence  between  the  scattering  angle  and 
the  impact  parameter  of  the  skewed  opacity  function  P(b)  in 
a  direct  collision.8,11  Such  behaviors,  a  rebound  mechanism 
near  the  threshold  and  shifting  to  peripheral  dynamics  at 
higher  Ec,  are  nicely  borne  out  in  recent  crossed-beam  stud¬ 
ies  on  C1+CH4  and  CD4,8,9  as  well  as  on  C1+CHD3.10  More¬ 
over,  closer  comparisons  of  the  H-  and  D-atom  abstraction 
channels  revealed  a  striking  isotope  effect.  At  Ec 
~7  kcal/mol,  for  example,  while  the  angular  distribution  of 
C1  +  CD4  reaction  peaks  near  90°  with  an  intensity  ratio  of 
/(180°)//(90°)  ~0.7,  that  of  C1+CH4  is  so  sharply  peaked  at 
60°  that  more  than  90%  of  product  fluxes  are  confined  within 
the  angular  range  of  30° -90°  with  a  ratio  of 
/( 1 80°)  //(60°)  ~  0.07  !8  The  Born-Oppenheimer  potential 
energy  surfaces  of  the  two  isotopic  reactions  are  identical; 
yet,  their  product  angular  distributions  differ  so  dramatically. 
Both  the  zero-point  energy  and  the  effect  of  tunneling 
through  the  centrifugal- shifted  barrier  to  reaction  are  in  favor 
of  H-atom  abstraction  channel,  which  have  been  proposed  to 
qualitatively  account  for  the  striking  isotope  difference  in 
angular  distributions  by  enhancing  the  characters  of  periph¬ 
eral  dynamics  when  a  H  atom  is  abstracted.8  Replacing  the 
tetrahedral  CH4/CD4  reactants  by  CH2D2  breaks  the  reaction 
symmetry,  as  such  the  Cl  atom  will  feel  different  peripheral 
interactions  as  the  H  atom  or  D  atom  is  attacked,  offering  a 
more  stringent  test  of  the  peripheral  mechanism. 

Also  closely  related  to  the  present  study  are  the  two 
remarkable  reports  on  the  bond  and  mode  selectivity  of  the 
title  reaction  with  vibrationally  excited  CH2D2  by  Zare  and 
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co-workers. 12,13  They  found  that  excitation  of  the  first  C-H 
overtone  of  CH2D2  promotes  only  the  H-atom  abstraction 
reaction,  and  excitation  of  the  first  C-D  overtone  of  CH2D2 
leads  preferentially  to  the  C-D  bond  cleavage.  In  addition, 
excitation  of  different  modes,  either  two  quanta  in  one  C-H 
oscillator  or  a  local  mode  containing  one  quantum  each  in 
two  C-H  oscillators,  produces  vastly  different  product  states 
and  angular  distributions.  The  major  findings  of  these  bond 
and  mode- selective  investigations  can  be  summarized  by  an 
intuitively  simple  spectator  model,12,13  in  which  the  Cl  atom 
interacts  with  the  single  C-H  oscillator  and  the  remainder  of 
the  methane  molecule  does  not  participate  in  the  reaction. 

To  complement  the  previous  studies  and  to  shed  more 
light  into  the  fascinating  mode  and  bond  selectivity  in  gen¬ 
eral,  we  report  here  an  extensive  study  on  reactions  of  Cl 
+  CH2D2^CHD2+HC1  and  CH2D  +  DC1.  Both  the  integral 
and  differential  cross  sections  were  measured  over  a  wide 
range  of  collision  energies,  —2-22  kcal/mol.  In  addition, 
the  reactivity  of  the  spin-orbit  excited  Cl*(2P1/2)  toward 
CH2D2(u  =  0)  and  that  of  bend/torsion-excited  CH2D2  with 
Cl (2P1/3)  are  characterized.  This  rather  extensive  data  set 
then  lays  a  firm  foundation  for  understanding  an  on-going 
study  in  this  laboratory  on  the  mode-selective  reactivity  with 
one  quantum  excitation  of  two  different  CH2- stretching 
modes  (symmetric  Vi  =  l  or  antisymmetric  v6=l)  in  Cl 

+ch2d2.14 

II.  EXPERIMENT 

The  experiments  were  carried  out  using  the  crossed- 
beam  apparatus  described  in  detail  elsewhere.8-10,15  In  brief, 
a  pulsed  high-voltage  discharge  method  was  used  to  generate 
the  Cl-atom  beam  (3.5%  Cl2  seeded  in  helium  at  6  atm  of 
total  pressure).  A  diluted  CH2D2  (—18%  CH2D2  seeded  in 
H2  at  6  atm)  from  an  E-L  valve16  was  used  to  cover  the 
entire  £c-range  of  this  study,  —2  to  greater  than  20  kcal/mol. 
The  reaction  products,  CHD2  and  CH2D,  were  probed  by 
(2+1)  resonance-enhanced  multiphoton  ionization  (REMPI) 
process17,18  using  a  time-sliced  ion  velocity-imaging  tech¬ 
nique  to  interrogate  the  product-pair-correlated 
information.15,19  REMPI  spectra  of  both  reaction  products 
indicated  predominant  formation  of  the  vibrational  ground 
state;  thus,  this  report  will  focus  on  the  pair-correlated  dy¬ 
namics  with  respect  to  the  CHD2(l>  =  0)  and  CH2D(u  =  0) 
products.  Experiments  on  probing  the  vibrationally  excited 
methyl  products  are  in  progress;  the  results  will  be  reported 
in  the  future. 

At  a  given  collision  energy,  the  time- sliced  image  was 
acquired  with  the  probe  laser  frequency  fixed  at  the  peak  of 
the  0°  Q-head  of  the  X2Bx-+ 3p2Bx  transition,  namely,  at  the 
two-photon  frequencies  (in  vacuum)  of  59  921  cm-1  for 
CHD2  and  of  59  940  cm'1  for  CH2D.17’18  The  observed 
g-head  is  substantially  broader  than  the  probe  laser  band¬ 
width.  Previous  investigations  from  this  laboratory  showed 
some  subtle  differences  in  pair-correlated  information  from 
two  alternative  modes  of  operation,  that  is,  either  by  fixing 
the  probe  laser  near  the  peak  of  the  2-head  or  by  scanning 
the  laser  frequency  back  and  forth  across  the  2-head  while 
acquiring  the  image.  ’  The  observed  differences  have  also 
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FIG.  1.  (Color)  Representative  time-sliced  raw  images  of  CHD2(00)  prod¬ 
ucts  at  four  different  collision  energies.  The  state  labelings  (i/  =  0),  (i/  =  l), 
(i/=0)^,  and  (i/  =  l)*  correspond  to  the  correlated  HC1  vibration  states  from 
the  reactions  of  (1)  C1+CH2D2(00)^CHD2(00)+HC1(u=0),  (2)  Cl 
+  CH2D2(00)  — ►  CHD2(00) + HCl(u  =  1 ) ,  (3)  Cl+CH2D2(ubend=l) 

^CHD2(00)+HC1(u  =  0),  and  (4)  Cr(2P1/2)  +  CH2D2(00)  -^CHD2(00) 

+HCl(i/  =  1),  respectively. 


been  demonstrated  to  originate  from  a  surprisingly  sensitive 
dependence  of  the  vibrational  branching  ratio  of  the  coprod¬ 
ucts  [HC1  (Ref.  21)  or  HF  (Ref.  20)]  on  the  probed  rotational 
A- states  of  the  methyl  products.  Although  the  vibrationally 
correlated  information  arguably  could  be  less  biased  toward 
the  probed  A- states  by  scanning  the  laser  frequency,  this 
mode  of  operation  generally  blurs  the  velocity-mapped  im¬ 
age  resolution.  Because  we  were  also  interested  in  the  reac- 
tivityofCl*(2P1/2)  and  of  the  bend-excited  CH2D2,  for  which 
the  velocity  resolution  is  of  prime  concern,  we  chose  to  fix 
the  laser  frequency  near  the  peak  of  2-head.  Thus,  the  pair- 

correlated  information  in  this  report  refers  to  that  with  re- 

20  21 

spect  to  the  v=0,  low  A  states  of  the  methyl  products.  ’ 

Depending  on  the  signal  strength,  each  time- sliced  im¬ 
age  was  taken  at  20  Hz  for  1-4  h,  accumulating  1  X  105  to 
6  X  105  total  events.  To  normalize  the  images  at  different 
collision  energies  for  obtaining  the  reactive  excitation  func- 
tions,  the  procedure  described  elsewhere  was  followed  in 
separate  measurements.  Relative  detection  sensitivities  of 
probing  the  0q  bands  of  CHD2  and  CH2D  are  yet  to  be  cali¬ 
brated.  Hence,  the  comparison  of  the  relative  cross  sections 
of  the  two  isotope  channels  was  merely  based  on  the  relative 
signal  strengths  under  otherwise  identical  experimental  con¬ 
ditions. 

III.  RAW  IMAGES  AND  PAIR-CORRELATED  ANGULAR 
DISTRIBUTIONS 

Figure  1  presents  typical  raw  images  of  the  CHD2(u 
=  0)  product  at  four  different  collision  energies,  and  simi¬ 
larly,  Fig.  2  for  the  CH2D(i/  =  0)  product.  The  energetics  of 
the  reaction  are  well  established:  AHrx=1.5  and  1.9  kcal/mol 
for  the  HC1+CHD2  and  DC1+CH2D  isotope  channels,  re¬ 
spectively;  so  as  the  vibrational  energies  of  HC1  (8.24  kcal/ 
mol  for  v  =  1)  and  DC1  (5.79  and  11.79  kcal/mol  for  v  =  1  and 
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FIG.  2.  (Color)  Similar  to  Fig.  1,  but  for  the  CH2D(00)+DC1  product 
channel. 


0  45  90  135  0  45  90  135  180 


ecm  (degree) 

FIG.  3.  The  pair-correlated  angular  distributions  of  the  C1+CH2D2 
-^CHD2+HC1  reactions  derived  from  the  four  images  shown  in  Fig.  1.  The 
line  symbols  are  the  same  for  all  four  panels  with  the  state  labels  (right- 
bottom  panel)  as  Fig.  1.  For  clarity,  the  intensities  of  the  minor  channels  are 
amplified. 


2,  respectively)  and  the  spin-orbit  energy  of  2.52  kcal/mol 
for  Cl*(2/>1/2).  Hence,  on  energetic  grounds,  all  ringlike  fea¬ 
tures  on  the  image  can  readily  be  identified  and  assigned  to 
the  vibrational  states  of  the  coproducts  from  different  reac¬ 
tant  states  as  labeled  in  the  figure  and  explained  in  the  cap¬ 
tion.  The  red  line  in  each  image  indicates  the  relative  veloc¬ 
ity  vector  in  the  center-of-mass  (c.m.)  frame  with  the  arrow 
pointing  to  the  0°,  which  is  defined  as  the  direction  of  the 
CH2D2  beam  in  the  c.m.  frame.  As  can  be  seen,  the  angular 
distribution  of  the  CHD2  or  CH2D  radical  coincidently 
formed  with  HCl(u  =  0)  or  DCl(u  =  0),  which  is  the  most 
prominent  feature  in  Fig.  1  or  Fig.  2,  evolves  from  backward 
near  threshold  to  sideways  and  then  into  forward  hemisphere 
at  successively  higher  Ec.  Such  a  systematic  angular  evolu¬ 
tion  with  Ec  appears  at  a  faster  pace  for  the  CHD2+HC1 
product  channel  than  for  CH2D  +  DC1. 

The  raw  images  shown  in  Figs.  1  and  2  do  not  possess 
the  cylindrical  symmetry  around  the  relative  velocity  axis 
because  of  the  nonuniform  detection  dependency  on  the 
product  laboratory  velocity.  ’  After  correcting  for  the 
density-to-flux  transformation,15  quantitative  angular  distri¬ 
bution  of  each  distinct  ring  can  then  be  deduced.  The  results 
of  such  image  analysis  are  presented  in  Figs.  3  and  4  for 
Figs.  1  and  2,  respectively.  A  few  remarks  are  in  order.  First, 
for  the  dominant  product  states  labeled  (u  =  0),  the  general 
trends  of  the  angular  distributions  of  the  two  isotopic  product 
channels  are  alike  except  the  different  rates  of  angular  evo¬ 
lutions  on  Ec,  as  alluded  to  above.  A  closer  comparison  re¬ 
veals  some  subtle,  perhaps  significant,  differences  between 
the  two  product  channels.  At  £c~21.7  kcal/mol,  the  HC1 
+  CHD2  channel  indicates  finite,  albeit  small,  fluxes  near  0°, 
which  is  evident  from  the  ring  closure  in  Fig.  1,  whereas  the 
flux  along  the  forward  direction  is  negligibly  small  in  the 
DC1  +  CH2D  channel.  In  addition,  the  ratio  of  the  peak  inten¬ 
sity  to  the  flux  at  180°  is  clearly  larger  for  HC1+CHD2  than 
that  for  DC1+CH2D.  Similar  behavior  is  seen  from  the  com¬ 
parison  of  HC1+CHD2  at  7.6  kcal/mol  to  DC1+CH2D  at  11.7 
kcal/mol,  even  though  the  former  is  at  lower  Ec.  The  dispar¬ 


ity  in  the  sharpness  of  angular  peaks  of  the  two  isotope  chan¬ 
nels  is  reminiscent  of  the  previous  findings  of  the  H- versus 
D-atom  abstractions  from  the  C1+CH4  and  CD4  reactions.8 

Second,  a  few  minor  image  features  are  worth  noting. 
The  angular  distributions  of  product  ground- state  pairs  from 
the  bend-excited  reaction  (see  Sec.  IV)  resemble  the  respec¬ 
tive  ground-state  reactions  (see  Figs.  3  and  4).  It  implies  that 
the  bend-excited  reaction  follows  nearly  the  same  pathway  as 
the  ground-state  reaction.  Similar  results  were  found  previ- 
ously  in  reactions  of  other  isotopomers.  ’  Since  the  bend- 
excited  reactants  correlate  adiabatically  to  the  bend-excited 
methyl  products,25,26  vibrationally  nonadiabatic  transitions 
must  occur  at  the  entrance  channel  in  order  for  bend- state 
reactions  to  proceed  over  the  ground- state  reaction  pathway 
by  funneling  the  bending  energies  into  the  product  transla¬ 
tional  and  rotational  energies  (i.e.,  acting  as  a  transitional 
mode).26  The  same  interpretation  was  proposed  previously 
for  the  other  isotopically  analogous  reactions.25,26 


0cm  (degree) 


FIG.  4.  Similar  to  Fig.  3,  the  pair-correlated  angular  distributions  of  the 
C1+CH2D2^CH2D+DC1  reactions  derived  from  the  four  images  shown  in 
Fig.  2. 
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As  to  the  vibrationally  excited  HCl(u  =  l)  from  ground- 
state  reactants  (the  thin  solid  lines  in  Fig.  3),  the  correlated 
angular  distributions  are  distinctly  different,  featuring  a  clear 
forward  and  a  weak  yet  discernible  backward  peaks  super¬ 
imposed  on  a  nearly  isotopic  component — in  analogy  to  the 
other  isotopomers.9,25  In  general,  such  an  angular  distribution 
is  suggestive  of  some  time-delay  reaction  mechanisms,27-31 
and  most  likely,  in  the  present  case,  a  short-lived  resonant 
complex  pathway.9,25,26  The  formation  of  DCl(u  =  1)  from  the 
ground- state  reactants,  however,  exhibits  a  rather  different 
angular  distribution  (Fig.  4),  which  appears  to  constitute  two 
components.  At  low  Ec,  the  distribution  resembles  the  corre¬ 
sponding  ground- state  product  pair  except  a  persistent  for¬ 
ward  component;  at  high  Ec,  it  turns  into  a  forward- 
backward  peaking  distribution  in  analogy  to  the  HCl(u  =  l) 
product.  In  other  words,  two  distinct  mechanisms  of  the  di¬ 
rect  and  the  resonant  reaction  pathways  dominate  over  the 
low  and  high  Ec  regimes,  respectively,  in  forming  the 
DC1(u  =  1)  +  CH2D(00)  product  pair.  The  angular  distribution 
for  the  DCl(u  =  2)  coproduct,  which  is  accessible  only  at  Ec 
>13.7  kcal/mol,  is  similar  to  DCl(u  =  l)  at  high  Ec,  i.e.,  a 
forward-backward  peaking  one. 

The  most  intriguing  finding  for  the  reaction  with  spin- 
orbit  excited  C r(2F1/2)  is  the  clear  formation  of  the  HC1/DC1 
coproducts  with  one  quantum  of  vibrational  excitation.  The 
correlated  angular  distribution  of  HCl(u  =  l)cl*  is  character¬ 
ized  by  the  distinct  forward-backward  asymmetric  peaks  (the 
dotted  lines  in  Fig.  3),  whereas  that  of  DC1(u  =  1)q*  exhibits 
additional  sideways  feature  that  resembles  the  ground- state 
pair  from  the  peripheral  pathway  (Fig.  4).  In  other  words,  the 
observed  angular  distribution  suggests  that  the  Cl*(2P1/2) 
+  CH2D2(u  =  0)  — >CHD2(00)+HC1(i;  =  1)  reaction  proceeds 
mainly  via  the  resonant  complex  pathway.  However,  both  the 
resonant  complex  and  peripheral  pathways  contribute  to  the 
formation  of  CH2D(00)  +  DC1(u  =  1). 

IV.  REACTIVITY  OF  BEND-EXCITED  CH2D2 

The  above  assignment  of  a  ringlike  feature  to  the 
(00,0)^  product  pair  from  the  reaction  with  bend-excited 
CH2D2  needs  further  clarification.  The  reactant  CH2D2  has 
nine  vibrational  degrees  of  freedom;  five  of  them  are  low- 
frequency  modes  (i.e.,  <1500  cm-1):  CD2  scissor 

(1033  cm-1),  CH2  rock  (1090  cm-1),  CH2  wag 
(1234  cm-1),  CH2  twist  (1333  cm-1),  and  CH2  scissor 
(1436  cm-1).32  The  first  two  modes  contain  merely 
—0.5  kcal/mol  more  energy  than  the  spin-orbit  excited 
Cl(2P1/2)  atom  (852  cm-1).  The  latter  must  also  be  present 
in  our  discharge-generated  beam,  as  evidenced  by  the 
(00, 1)C1*  feature  as  well  by  the  previous  finding  in  C1  +  CH4.9 
Therefore,  energetic  argument  alone  is  not  sufficient  for  an 
unambiguous  assignment.  The  unequivocal  evidence  came 
from  the  separate  experiment  of  using  a  heated  CH2D2 
pulsed  valve,  as  first  demonstrated  in  the  studies  of  Cl 
+CH4/CD4  (Ref.  24)  and  0(3P)+CD4/CHD3  (Ref.  33)  re- 
actions. 

Shown  in  Fig.  5  are  the  results  from  a  back-to-back  ex¬ 
periment  of  the  HC1(u)  +  CHD2(00)  channel,  for  which  the 
CH2D2-valve  temperature  was  varied.  The  molecular-beam 


FIG.  5.  (Color)  Effect  of  low-lying  hot  bands  of  CH2D2  reactants  on  the 
observed  images.  The  upper  panels  display  the  CHD2(00)  product  images  at 
two  different  CH2D2  pulsed  valve  temperatures.  The  molecular-beam  inter¬ 
section  angles  were  adjusted  slightly  so  that  the  two  images  were  acquired  at 
the  same  £c~  5.8  kcal/mol.  The  analysis  of  product  speed  distributions  are 
presented  in  the  lower  panel.  The  two  vertical  arrows  mark  the  maximal 
speeds  from  reactions  with  the  ground-state  and  bend-excited  reactants,  re¬ 
spectively.  For  ready  comparison,  the  intensities  of  two  distributions  are 
normalized  by  the  peaks  from  the  ground-state  reaction  (the  slow  compo¬ 
nent);  significantly  enhanced  faster  component  at  higher  temperature  is  viv¬ 
idly  displayed. 

intersection  angles  were  adjusted  slightly  so  that  the  colli- 
sional  energies  from  different  source  temperatures  remained 
about  the  same,  Ec~  5.8  kcal/mol.  Even  a  casual  inspection 
of  the  raw  image  reveals  a  significant  enhancement  in  the 
(u  =  0)*  feature  at  higher  temperature  (439  K).  More  quanti¬ 
tative  comparison  is  demonstrated  from  the  fully  analyzed 
product  speed  distribution  P(u),15  which  shows  two  well- 
resolved  peaks  (the  lower  panel  of  Fig.  5).  For  ready  com¬ 
parison,  the  P(u)  distributions  from  the  two  source  tempera¬ 
tures  are  normalized  by  the  peak  height  of  the  slow 
component,  which  is  ascribed  to  C1+CH2D2(l’  =  0)^HC1(u 
=  0)  +  CHD2(00).  As  is  seen,  the  P{u)  distribution  for  439  K 
displays  a  more  intense  fast  component  with  the  intensity 
ratio  of  the  fast-to-slow  components  increasing  from  —9%  at 
320  K  to  —32%  at  439  K.  Because  the  Cl  beam  was  kept  the 
same,  this  difference  by  heating  the  CH2D2  valve  provided  a 
compelling  evidence  for  the  assignment  of  the  fast  compo¬ 
nent  mainly  to  the  (0,00)^  product  pair  from  the  reaction 
with  thermally  excited  CH2D2. 

Previous  studies  on  C1+CH4/CD4  reactions  demon¬ 
strated  that  although  supersonic  expansion  yielded  a  transla- 
tionally  and  rotationally  cold  (—5-10  K)  methane  beam, 
the  low-lying  vibration  manifolds  remain,  by  and  large,  ther- 
mally  populated.  Assuming  that  the  similar  situation  holds 
for  the  present  CH2D2  beam,  one  can  then  estimate  the  reac- 
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tivity  of  bend-excited  reactants.  The  above  fast-to-slow  ratio 
corresponds  to  /n0a0  with  nt  denoting  the  population 

of  the  vibrationally  excited  state  i  of  CH2D2  in  the  beam  and 
of  being  the  corresponding  reaction  cross  section  (the  sub¬ 
script  “0”  indicates  the  vibrational  ground- state  reactant). 
Since  the  thermally  populated  vibration  manifolds  follow  the 
Boltzmann  distribution  (i.e.,  without  state  selection)  and,  en¬ 
ergetically,  the  P(u)  distribution  of  the  fast  component  rules 
out  the  contributions  from  higher-lying  vibration  levels, 
therefore  only  the  five  low-frequency  modes  need  to  be  con¬ 
sidered.  We  further  approximated  /n0o0  as 

o-#E^-/n0(r0,  with  o*  representing  the  average  cross  section 
over  those  five  modes.  The  Boltzmann  factors  at  320 

and  439  K  are  0.025  and  0.102,  respectively.  Hence,  the 
estimated  vibrational  enhancement  factors  <x#/ob  are  3.6  (or 
0.09/0.025)  and  3.1  (or  0.32/0.102)  from  the  analysis  of  the 
data  at  two  temperatures.  It  is  gratifying  to  note  the  self- 
consistency  of  the  two  estimated  enhancement  factors,  which 
strongly  suggests  that,  within  our  experimental  uncertainties, 
the  contribution  of  Cl*  +  CH2D2(u  =  0)  to  the  formation  of 
HC1(u  =  0)  +  CHD2(00)  is  relatively  minor  at  this  collision  en¬ 
ergy. 

Similarly,  Fig.  6  presents  the  images  and  the  P(u)  distri¬ 
butions  for  the  other  isotope  channel,  DC1+CH2D.  The  ratios 
of  fast-to-slow  components  are  —8%  and  —30%  at  320  and 
436  K,  respectively,  which,  in  turn,  yield  the  respective 
<x#/<x0  — 3.2  and  —3.0  for  the  bend-excited  enhancement 
factor  at  Ec  — 5.8  keal/mol;  again,  the  agreement  was  quite 
reasonable.  It  is  interesting  to  note  that  the  vibrational  en¬ 
hancement  factor  of  —3  for  bend-excited  CH2D2  in  both 
isotope  channels  is  virtually  the  same  as  the  factors  obtained 
for  the  bend-excited  C1  +  CH4/CD4  (Refs.  24  and  34)  and 
CHD3  (Ref.  25)  reactions. 
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FIG.  7.  The  relative  reactivity  of  the  four  state-to- state  reactions  listed  in  the 
caption  of  Fig.  1  when  probing  the  CHD2(00)  products.  The  left-top  panel 
shows  the  relative  excitation  functions  weighed  by  the  respective  reactant- 
state  concentrations.  The  remaining  three  panels  present  the  reactivity  ratios 
to  the  dominant  ground-state  reaction  cross  section  cr0.  To  obtain  the  relative 
cross  section  of  bend-excited  reactants,  a^/(T0,  the  ordinate  of  the  left- 
bottom  panel  should  be  divided  by  n*/n0,  which  is  about  0.025  at  the 
CH2D2  source  temperature  of  320  K.  The  vertical  arrows  in  two  right  panels 
indicate  the  respective  energetic  thresholds  of  the  two  reaction  channels. 


V.  CORRELATED  EXCITATION  FUNCTIONS 
AND  RELATIVE  REACTIVITY 

After  normalizing  each  individual  methyl  product  image 
taken  at  different  collision  energies  and  correcting  for  the 
density-to-flux  transformation,  the  energy  dependence  of  re¬ 
active  fluxes  on  the  respective  product-pair  channels  can 
then  be  deduced.22  The  resultant  state-to- state  excitation 
functions  for  the  HC1  +  CHD2  channel  are  presented  in  Fig.  7. 
As  discussed  above,  four  pair-correlated  reaction  channels 
were  identified  from  the  product  image  of  CHD2(00).  In  ad¬ 
dition  to  the  two  product  pairs  CHD2(00)  +  HC1(u  =  0  and  1) 
from  the  ground- state  reaction,  the  contributions  from  the 
excited  reactants,  i.e.,  Cl(2P3/2)  +CH2D^(bend) 

^CHD2(00)  +  HC1(u=0)  and  C1*(2P1/2)+CH2D2(i;=0) 
^CHD2(00)  +  HC1(u  =  1),  were  observed  and  analyzed.  The 
relative  signal  of  the  latter  two  obviously  depends  on  the 
initial  population  (n)  of  the  excited  reactant  state  in  the  mo¬ 
lecular  beam.  Hence,  the  excitation  functions  shown  in  Fig.  7 
(left- top  panel)  are  plotted  as  ncr.  The  other  three  panels 
display  the  respective  ^-dependences  of  the  ratios  of  those 
minor  image  features  on  n0o0  of  the  dominant  reaction  chan¬ 
nel  Cl + CH2D2(u  =  0)  CH2D(00)  +  HCl(u = 0) . 

Compared  to  the  ground-state  reaction,  the  (00,0)^ 
product  pair  from  the  bend-excited  reactants  (Fig.  7,  left- 
bottom  panel)  displays  a  monotonic  decline  with  increasing 
Ec.  The  estimated  n^/n0  is  about  0.025  at  CH2D2  source 
temperature  of  320  K;  thus,  the  ordinate  for  o^  /  o0  should  be 
multiplied  by  40,  yielding  a  vibrational  enhancement  of  at 
least  a  factor  of  2.  It  is  tempting  to  attribute  this  enhance¬ 
ment  factor  to  the  energetic  effect  in  that  whereas  the 
ground-state  reaction  possesses  a  barrier  to  reaction,  the 
bend-excited  reactants  contain,  on  average,  additional 
—  3.5  keal/mol  of  energy  for  reaction.  If  the  initial  bend 
excitation  of  CH2D2  is  readily  coupled  to  the  reaction  coor¬ 
dinate  through  curvature  couplings  before  the  transition 
state,  the  transfer  of  bending  energy  will  then  help  surmount 
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the  reaction  barrier,  promoting  the  reaction  rate.  Yet,  the  ob¬ 
served  vibrational  enhancement  becomes  more  pronounced 
at  lower  Ec ,  which  is  in  sharp  contrast  to  the  previous  finding 
of  a  nearly  ^-independent  enhancement  (at  low  Ec)  for 
bend-excited  CHD3  in  the  C1+CHD3— >CD3(00)  +  HC1 
reaction.25  The  observation  of  different  /^-dependencies  of 
bend-mode  enhancements  of  two  isotope  reactants  is  puz¬ 
zling,  perhaps  points  to  the  subtle  differences  in  the  curva¬ 
ture  couplings  and/or  the  Coriolis  interactions  among  the 
bending/torsional  modes.  Further  investigation  is  warranted. 

The  concomitant  HC1  vibrational  branching  ratio  crl/  cr0 
from  the  ground-state  reaction  (Fig.  7,  right-top  panel,  the 
vertical  arrow  indicates  the  energetic  threshold)  remains 
small  even  at  higher  Ec ,  implying  that  the  reaction  of  ground- 
state  reactants  is,  by  and  large,  vibrationally  adiabatic.  Same 
conclusion  was  drawn  in  the  previous  studies  of  C1  +  CH4 
(Ref.  9)  and  CI+CHD3.25  As  to  the  CHD2(00)  +  HCl(n  =  1) 
pair  from  the  spin-orbit  excited  C I  ' ( 2  /J ,  /2 )  reaction  (Fig.  7, 
right-bottom  panel),  the  relative  reactivity  n*a\/n0a0  exhib¬ 
its  a  threshold  around  Ec=13  kcal/mol,  as  anticipated  on 
energetic  ground  (the  arrow).  The  relative  population  of 
Cl*(2F1/2)  in  the  beam,  n* ln0,  was  not  measured  in  this 
study,  but  is  expected  to  be  less  than  1/2  (i.e.,  the  relative 
degeneracy  factors  of  the  excited  2P1/2  t0  the  ground  2F3/2 
states);  otherwise,  population  inversion  ensues  from  dis¬ 
charging  Cl2  in  He.  With  this  caveat,  the  relative  reactivity  of 
ci*(2p1/2)  toward  CH2D2  appears  quite  significant  at  higher 
collision  energies,  e.g.,  a\l a0  will  approach  9%  at  Ec 
—  20  kcal/mol  or  likely  higher  if  n*/n0  is  significantly  less 
than  1/2.  Perhaps  even  more  remarkable  is  the  fact  that  the 
concomitant  product  state  here  refers  to  the  vibrationally  ex¬ 
cited  HCl(u  =  l);  no  image  feature  can  be  unambiguously 
ascribed  to  the  energetically  more  readily  accessible  HCl(u 
=  0)  coproduct,  as  pointed  out  in  Sec.  IV.  It  is  interesting  to 
note  that  the  formation  of  CP(2P1/2)  product  from  the  reverse 
reaction  of  CH3/CD3+HC1  at  Ec  —  20  kcal/mol  was  re¬ 
cently  reported. 35-37  However,  in  that  study  the  CH3/CD3 
reactant  was  generated  photolytically  from  CH 3//  CD37,  thus 
containing  a  significant  amount  of  internal  energy,  whereas 
the  HC1  reactant  is  mostly  in  the  v  =  0  state.  Because  of  the 
differences  between  the  initially  prepared  reactant  states  of 
that  study  and  the  probed  product  states  of  this  study,  it  is 
difficult  to  compare  the  two  investigations. 

Similar  results  for  the  other  isotope  channel  CH2D 
+  DC1  are  presented  in  Fig.  8.  For  the  bend-excited  reaction 
(left-bottom  panel),  again  the  relative  reactivity  n^a^  ln0a0 
declines  with  the  increase  in  Ec ,  and  the  estimated  <x#/c r0 
should  be  scaled  up  by  —40  to  account  for  the  small  relative 
population  of  n* ln0  in  the  CH2D2  beam.  Compared  to  the 
CHD2  +  HC1  channel,  the  vibrational  enhancements  from  the 
bend-excited  reactants  are  comparable  except  near  the 
ground-state  reaction  thresholds.  For  reaction  with  ground- 
state  reactants,  the  formation  of  DCl(u  =  2)  becomes  ener¬ 
getically  accessible  at  higher  Ec.  Indeed,  both  ax!  a -0  and 
a2/a0  become  finite  at  the  respective  energetic  thresholds 
(indicated  by  the  arrows).  Their  combined  reactivity,  how¬ 
ever,  remains  small  with  respect  to  the  ground-state  product 
pair.  Thus,  the  reaction  of  the  ground-state  reactants  for  this 
isotope  channel  is  also  largely  vibrationally  adiabatic,  al¬ 


■  0 

0  0  bend 
a  1  (X10) 
•  2  (X20) 

■ 

■ 

■  ■ 

■  A- 

▲ 

aAa 

▲ 

M.  A  A 

O  *  InnO 

O 

O 

0  0  0  0  0  o  0  0 

*• 

# 

# 

# 

5  10  15  20  5  10  15  20 

Ec  (kcal/mol ) 


0.06 

0.03 

0.00 

0.012 

0.006 

0.000 


o 


o 


b 


CN 


b 


FIG.  8.  Same  as  Fig.  7,  but  for  the  CH2D+DC1  product  channels. 


though  its  adiabaticity  does  not  appear  as  strong  as  the 
HC1+CHD2  channel.  Unfortunately,  due  to  the  beam  back¬ 
ground  (in  the  forward  direction)  problem  and  the  higher 
demand  in  image  resolution  for  this  isotope  channel,  only  for 
a  few  images  of  CH2D  over  a  limited  Ec  range  can  the  spin- 
orbit  excited  features  be  analyzed;  thus,  its  excitation  func¬ 
tion  cannot  be  reliably  determined. 

VI.  Ec-EVOLUTION  OF  PAIR-CORRELATED  ANGULAR 
DISTRIBUTIONS 

By  combining  the  pair-correlated  excitation  function 
presented  in  Sec.  V  and  the  angular  distribution  at  each  col- 
lisional  energy  as  illustrated  in  Sec.  Ill,  a  three-dimensional 
dal  d£l-0-Ec  plot  can  be  obtained.6,19  Such  a  representation 
has  proven  to  provide  a  global  and  illuminating  view  of  the 
underlying  reaction  mechanism,  in  particular  for  the  “discov¬ 
ery”  of  reactive  resonances  in  the  isotopically  analogous  re¬ 
actions  of  C1+CH4  (Ref.  9)  and  CHD3.25  The  results  of  the 
present  four  product  pairs  when  CHD2(00)  was  probed  are 
summarized  in  Fig.  9.  As  is  seen,  the  patterns  of  (uHC1=0) 
and  (uHci=0)^  are  alike,  displaying  a  slowly  evolving  ridge 
structure  that  runs  from  a  weak  backward  feature  at  low  Ec 
toward  an  intense  peak  in  the  forward  hemisphere  at  higher 
Ec.  Virtually  identical  patterns  were  noted  previously  for  the 
analogous  product  pairs  from  the  C1+CH4  (Ref.  9)  and 
CHD3  reactions,25  which  were  ascribed  to  a  manifestation  of 
a  direct  reaction  governed  by  rebound  mechanism  near 
threshold  and  then  shifting  to  peripheral  dynamics  at  higher 
E^8"10 

However,  the  two  patterns  shown  in  the  lower  panel  are 
distinctly  different.  For  (uHC1=l),  i.e.,  the  HCl(u  =  l) 
+  CHD2(00)  product  pair  from  the  ground-state  reactants, 
clear  forward-backward  peaking  distributions  are  seen  in  ad¬ 
dition  to  the  ridge  structure.  A  closer  inspection  further  un¬ 
veils  that  the  ridge  structure  here  does  not  look  like  the 
above  two  for  (uHC1=0)  and  (uHC1=0)^.  We  surmised  that  the 
ridge  here  actually  comprises  two  overlapped  ridges:  a  fast- 
evolving  (with  respect  to  collision  energy)  one  with  higher 
backward  intensity  at  lower  Ec,  and  as  Ec  increases,  it 
merges  into  a  slow-evolving  one  that  is  similar  to  the  ones  in 
the  upper  panel.  The  fast-evolving  ridge  near  the  energetic 
threshold,  together  with  the  forward-backward  peaking  dis- 
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FIG.  9.  (Color)  Three-dimensional 
representations  of  dal  d£l-d-Ec  that 
summarize  the  evolution  of  correlated 
angular  distributions  with  collisional 
energies  for  reactions  of  C1+CH2D2 
->CHD2(00)+HC1(i;).  Note  the  differ¬ 
ent  intensity  scales  of  each  panel.  The 
angular  distributions  are  weighed  by 
their  respective  populations  of  reactant 
states  (i.e.,  similar  to  the  left-top  panel 
of  Fig.  7).  Distinct  patterns  are  clearly 
seen,  suggestive  of  two  dominant  un¬ 
derlying  reaction  pathways:  a  direct 
pathway  evolving  from  the  rebound 
mechanism  at  low  Ec  to  a  peripheral 
dynamics  at  high  Ec  for  (uHci=0)  and 
(uHC]  =  0)^,  and  an  additional  resonant 
complex  mechanism  also  contributes, 
leading  to  the  product  channels  la¬ 
beled  as  (uHC1=l)  and  (uHC1=l)*. 


tribution  at  higher  Ec ,  is  reminiscent  of  the  characteristic 
pattern  for  reactive  resonance  uncovered  in  the  F+HD 
^HF+Z)  reaction.38-40  The  physical  origin  of  this  resonance 
signature,  namely,  how  a  dynamical  resonance  gives  rise  to  a 
fast-evolving  ridge  followed  by  forward-backward  peaking 
angular  distribution  in  bimolecular  reaction,  has  been  eluci¬ 
dated  and,  to  a  large  extent,  understood.6,39-41  Based  on  pat¬ 
tern  comparison,  we  assert  that  the  formation  of  (uHC1=l) 
from  the  ground-state  reactant  invokes  both  the  resonance 
and  the  direct  peripheral  reaction  pathways,  which  is  also 
what  we  concluded  previously  for  forming  HCl(u  =  l)  in  re¬ 
actions  of  C1  +  CH4  (Ref.  9)  and  CHD3.25  A  recent  quantum 
dynamics  and  quasiclassical  trajectory  study  of  C1+CH4  on  a 
pseudotriatomic  ab  initio  based  surface  also  unveiled  similar 
imprints  of  a  Feshbach  resonance  for  the  production  of 

HCl(u  =  l)  (Ref.  42)  in  support  of  our  experimental 

9 

suggestion. 

As  to  the  spin-orbit  excited  reaction  of  Cr(2JP1/2) 
+  CH2D2^HC1(u  =  1)  +  CHD2(00),  labeled  as  (uHci=l)*  in 
Fig.  9,  the  pattern  is  dominated  by  the  forward-backward 
peaking,  with  the  ridge  structures  (both  the  low-energy  reso¬ 
nance  ridge  and  the  high-energy  peripheral  ridge)  being 
nearly  invisible.  Adiabatically,  the  spin-orbit  excited 
Cl*(2F1/2)  cannot  react  with  methane  at  the  energy  range  of 
this  study.  Thus,  the  formation  of  (uHci-l)*  must  proceed 
through  an  electronically  nonadiabatic  process.  The  general 
pattern  shown  in  Fig.  9  is  reminiscent  of  that  for  F*(2P  1/2) 
+  HD— >HF(u  =  3)  +  D,43  which  has  been  ascribed  to  a 
resonance-mediated  (spin-orbit)  nonadiabatic  process.  Based 
on  the  prior  ab  initio  calculations  of  the  C1+  methane 
reactions44-46  and  the  vibrationally  adiabatic  correlation 
scheme,9,25,26  a  dynamical  well  that  correlates  the  reactant- 
pair  C1(2jP3/2)  +  CH2D2(l’1  =  1)  to  the  product-pair  HCl(u  =  l) 
+  CHD2(00)  is  anticipated.  We  speculate  that  the  nonadia¬ 
batic  deactivation  from  the  spin-orbit  excited  Cl*(2P1/2) 


+  CH2D2  to  the  ground  Cl(2P3/2)  +  CH2D2  surfaces  in  the  en¬ 
trance  channel,  which  is  a  prerequisite  for  Cr(2P1/2)  to  be 
reactive,  is  likely  accompanied  by  a  simultaneous  promotion 
of  the  CH2-symmetric  stretch  (rq  =  1)  of  the  CH2D2  reactant. 
Asymptotically,  one  quantum  excitation  of  the  rq-mode  re¬ 
quires  —  8.5  kcal/mol,  which  is  significantly  larger  than  the 
spin-orbit  energy  of  Cl*(2P1/2),  2.52  kcal/mol.  However,  the 
ab  initio  calculations  indicated  that  the  vibrational  frequency 
of  the  urmode  will  decrease  drastically  upon  the  approach 
of  the  Cl  atom  due  to  the  strong  curvature  coupling  to  the 
reaction  coordinate,44-46  which  would  narrow  the  energy  gap 
between  the  two  vibronic  surfaces  along  the  reaction  coordi¬ 
nate,  facilitating  the  suggested  (vibrationally  mediated)  elec¬ 
tronic  surface  hopping  process.  At  sufficiently  high  Ec ,  it 
could  then  lead  to  the  observed  HCl(u  =  1)  +  CHD2(00)  prod¬ 
uct  pair  through  the  reactive  resonance  supported  by  the  dy¬ 
namic  well  on  the  (rq  =  l)  vibrationally  adiabatic  curve. 

The  qualitative  interpretation  given  here  is  obviously 
tentative  and  does  not  exclude  other  lower  frequency  vibra¬ 
tional  modes  of  CH2D2  that  are  much  closer  in  resonance  to 
the  spin-orbit  splitting  of  the  Cl  atom  from  participating  in 
the  electronically  nonadiabatic  dynamics.  (In  fact,  as  men¬ 
tioned  above,  some  bending/torsional  motions  of  CH2D2  are 
Corialis  coupled  to  the  symmetric-stretching  mode.)  We 
merely  wish  to  emphasize  the  essential  role  of  the 
symmetric- stretching  mode  in  facilitating  the  spin-orbit 
nonadiabatic  transition  in  order  to  account  for  the  distinct 
formation  of  the  HCl(u  =  1)  +  CHD2(00)  pair.  More  theoreti¬ 
cal  work  is  invited. 

Similarly,  Fig.  10  presents  the  result  of  the  four  product 
pairs  observed  in  the  CH2D  +  DC1  isotopic  channel.  Although 
the  details  are  not  quite  the  same  as  the  CHD2+HC1 
channel — for  example,  the  formation  of  (uDC1=0)^  from  the 
bend-excited  reaction  shows  a  pronounced  direct  swath38,40 
in  the  backward  hemisphere,  the  gross  features  of  the  periph- 
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FIG.  10.  (Color)  Same  as  Fig.  9,  but  for  reactions  of 
Cl+CH2D2^CH2D(00)+DCl(u).  Due  to  the  beam 
background  and  image  resolution  problems,  the  results 
for  spin-orbit  excited  reactivity  cannot  be  obtained  reli¬ 
ably  over  a  wide  range  of  energy;  thus  are  not  shown. 
The  patterns  indicate  that  the  correlated  ground- state 
DC1  products,  (uDci=0)  and  (uDci=0)^  are  formed  pre¬ 
dominantly  by  direct  rebound/peripheral  reaction 
mechanisms,  whereas  the  reactive  resonance  pathway 
also  contributes  to  the  formation  of  vibrationally  ex¬ 
cited  DC1,  that  is,  (uDC1=l)  and  (uDC1=2). 


eral  ridge  (the  upper  panel)  and  the  resonance  characters  (the 
lower  panel)  are  retained.  Hence,  we  assert  that  the  forma¬ 
tions  of  (uDC1=0)  and  (uDC1=0)^  are  of  direct  reaction 
mechanism,  whereas  the  productions  of  (uDC1=l)  and  (uDC1 
=  2)  proceed  through  both  the  direct  and  resonance  pathways. 
Judging  from  the  appearance,  the  formation  of  DCl(u  =  2) 
+  CH2D(00)  seems  to  be  dominated  more  by  the  resonant 
reaction  mechanism. 

VII.  CORRELATED  ENERGY  DISPOSAL 

It  is  instructive  to  look  back  on  a  more  average  quantity, 
the  correlated  energy  disposal  of  the  two  isotopic  product 
channels  from  the  dominant  ground-state  reaction,  to  gain 
further  insights  into  their  reactivity.  The  results  for 
C\(2Pj/2) + CH2D2(i>  =  0)  — >  HCl(i;)  +  CHD2(00)  are  summa- 
rized  in  Fig.  11.  The  fractional  energy  disposals  are  depicted 
in  the  upper  panel  and  the  absolute  energy  (in  kcal/mol)  are 
presented  in  the  lower  panel.  Note  the  different  scales  on  the 
right  and  left  ordinates.  Most  of  the  available  energy  appears 
as  product  translational  energy;  moreover  as  Ec  increases,  the 


FIG.  11.  Correlated  energy  disposal  of  the  C1+CH2D2(00)^CHD2(00) 
+HCl(i;)  reaction.  The  fractional  energy  disposal  is  presented  in  the  upper 
panel,  and  the  energy  disposal  (in  kcal/mol)  in  the  lower  panel.  Note  the 
different  scales  of  the  right  and  left  ordinates.  Predictions  from  a  simple 
kinematics  model  are  given  by  the  gray  lines. 


fractional  kinetic  energy  release  fT  remains  nearly  constant, 
—  88%.  The  remaining  is  largely  deposited  into  HC1  rotation, 
fR.  Note  that  the  small  HC1  vibrational  fraction  fv  appears  at 
the  expense  of  HC1  rotational  degree  of  freedom. 

Also  shown  as  the  heavy  gray  line  in  the  lower  panel  is 
the  prediction  of  the  kinetic  energy  release  ET  from  the  con¬ 
sideration  of  kinematic  constraints.8  The  title  reaction  is  es¬ 
sentially  a  heavy -flight-heavy  system,  involving  a  H-atom 
transfer  from  the  CHD2  moiety  to  the  attacking  Cl  atom.  For 
such  a  mass  combination  of  an  endothermic  reaction,  a 
simple  kinematics  model  predicts  that  the  average  kinetic 
energy  release  can  be  expressed  as  ET=(EC-  A#rx)cos2  ft.41 
Here,  A Hrx  is  the  endothermicity  and  f3  is  the  skew  angle, 
cos 2  /3=mAmc/mABmBC,  for  the  A+BC^AB  +  C  reaction. 
The  excellent  agreement  between  the  experimental  data  and 
the  model  prediction  over  such  a  wide  range  of  Ec,  in  which 
no  adjustable  parameter  is  invoked,  provides  a  compelling 
evidence  for  the  dominant  kinematics  role  in  governing  this 
attribute,  as  well  the  spectator  nature  of  the  methyl  moiety  in 
the  ground-state  reaction.  This  conclusion  reinforces  our  pre¬ 
vious  findings  for  the  C1+CH4/CD4  (Ref  8)  and  C1+CHD3 
reactions.10 

Figure  12  presents  similar  analysis  for  the  DC1+CH2D 


FIG.  12.  Same  as  Fig.  11,  but  for  the  Cl+CH2D2(Oo)^CH2D(00) 
+DCl(u)  reaction. 
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isotope  channel.  Once  again,  the  simple  kinematics  model 
accounts  for  the  results  remarkably  well.  A  closer  compari¬ 
son  between  the  two  isotope  channels  reveals  that  the  aver¬ 
age  kinetic  energy  release  ET  for  the  H-atom  abstraction 
channel  is  slightly  larger  than  that  for  the  D-atom  channel. 
Alternatively,  this  may  be  more  apparent  from  the  compari¬ 
son  of  the  HC1/DC1  rotational  energy  release,  ER.  As  dis¬ 
cussed  above,  the  reaction  is  mainly  vibrationally  adiabatic 
in  the  sense  that  the  vibrational  ground-state  reactants  yield 
mostly  the  vibrational  ground-state  products.  Hence,  the  re¬ 
maining  energy  release  from  the  kinematics  model,  (Ec 

-  AF/rx)sin2  /3,  is  mainly  deposited  into  the  HC1/DC1  rotation. 
Indeed,  this  is  what  we  observed,  as  shown  by  the  thin  gray 
lines  in  Figs.  11  and  12.  For  the  two  isotopic  channels, 
£^(DC1)  is  about  twice  larger  than  ^(HCl);  for  example,  at 
Ec= 20  kcal/mol  ^(HCl)  and  £^(DC1)  are  1.5  and  2.8  kcal / 
mol,  respectively.  For  a  A+BC  reaction,  the  angular  momen¬ 
tum  disposal  can  be  expressed  as  j'  =1  sin2  (3+j  cos2  (3 
+ d  cos2  /3.47  Here,  j'(j)  is  the  product  (reactant)  rotational 
angular  momentum  (we  neglected  the  low-ZV  angular  mo¬ 
menta  of  the  probed  methyl  products,  which  were  kept  the 
same  for  all  Ec  in  this  study),  /  is  the  initial  orbital  angular 
momentum  of  the  colliding  reactants,  and  d  is  a  transforma¬ 
tion  vector  incorporating  the  effect  of  reaction  dynamics.47 
Since  j  ~  0  for  a  supersonic  CH2D2  beam  and,  as  just  shown, 
the  kinematics  dominates  the  energy  disposal,  we  conjec¬ 
tured  that  the  dynamical  factor  d  is  also  small.  In  other 
words,  we  approximated  j'  —/  sin2  /3= 0.081/  (or  0.16/)  for 
the  HC1  (or  DC1)  channel.  The  rotational  constants  for  HC1 
and  DC1  are  10.44  and  5.4  cm-1,  respectively.  From  the  av¬ 
erage  rotational  energies  ER  given  above,  one  has  j'HCi~ 7ft 
and  7qC1~13^,  which,  in  turn,  lead  to  /  —  86^  and  81ft  for 
the  HC1  and  DC1  channels,  respectively,  at  Ec 

—  20  kcal/mol. 

The  estimated  /-values  here  refer  to  the  average  orbital 
angular  momentum  of  the  two  reactants.  Because  the  reac¬ 
tion  is  of  the  peripheral  type,  the  opacity  function  strongly 
favors  larger  /.  Hence,  the  average  value  of  /  should  be  quite 
close  to  the  peak-/  of  the  opacity  function,  which  is  in  rea¬ 
sonable  accord  with  theoretical  results  for  the  isotopically 
analogous  reactions.3,10,48,49  Moreover,  the  closeness  of  the 
two  /-values  estimated  for  the  HC1  and  DC1  channels  implies 
that  the  respective  reactive  cross  sections  cannot  be  very  dif¬ 
ferent.  From  Figs.  7  and  8  (left-top  panels),  the  measured 
cross  sections,  without  accounting  for  the  detection  sensitiv¬ 
ity  differences  when  probing  the  low  A-states  of  CHD2(0q) 
and  CH2D(0°),  are  nearly  the  same.  It  then  suggests  that  the 
(2+1)  REMPI  detections  of  the  two  isotopic  methyl  radicals 
may  have  approximately  the  same  sensitivity — a  point  needs 
future  confirmation. 

VIII.  CONCLUSIONS 

In  this  report,  we  focus  on  the  correlated  reaction  dy¬ 
namics  associated  with  the  ground  state  of  methyl  radicals 
CHD2(00)  and  CH2D(00),  which  are  the  major  products  of 
the  title  reaction.  Thanks  to  the  high  resolution  of  the  time- 
sliced  velocity-imaging  detection  scheme,  the  raw  images 
revealed  the  coincidently  formed  HC1/DC1  coproducts  in  a 


state-resolved  manner  from  not  only  the  ground-state  reac¬ 
tants  but  also  the  bend-excited  CH2D2  and  the  spin-orbit  ex¬ 
cited  a*(2P1/2)  reactants.  In  other  words,  state-to-state  cor¬ 
related  dynamics  for  three  different  combinations  of  reactant 
pairs,  i.e.,  Cl(2P3/2)  +  CH2D2(00),  Cl(2F3/2)  +  CH2D2(L>bend 
=  1)#,  and  Cl*(2P1/2)  +  CH2D2(00),  were  acquired  simulta¬ 
neously  by  imaging  the  CHD2(00)  and  CH2D(00)  products. 
The  relative  reactivity  of  the  observed  eight  different 
product-pair  channels  were  characterized,  where  possible, 
and  their  dynamical  attributes  were  elucidated.  Two  distinct 
reaction  mechanisms,  direct  peripheral  and  reactive  reso¬ 
nance,  were  put  forward  to  account  for  the  observations. 
Similar  isotope  effects  as  those  found  previously  in  Cl 
+  CH4/CD4  and  C1+CHD3  were  also  observed,  which  appear 
to  be  a  universal  trait  of  this  benchmark  reaction  of  chlorine 
atom + methane. 
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